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ABSTRACT 

Context. The solar photospheric oxygen abundance is still widely debated. Adopting the solar chemical composition based on the 
‘low’ oxygen abundance, as determined with the use of three-dimensional (3D) hydrodynamical model atmospheres, results in a well- 
known mismatch between theoretical solar models and helioseismic measurements that is so far unresolved. 

Aims. We carry out an independent redetermination of the solar oxygen abundance by investigating the center-to-limb variation of the 
O i IR triplet lines at 777 nm in different sets of spectra. 

Methods. The high-resolution and high signal-to-noise solar center-to-limb spectra are analyzed with the help of detailed synthetic 
line profiles based on 3D hydrodynamical C05B0LD model atmospheres and 3D non-LTE line formation calculations with NLTE3D. 
The idea is to exploit the information contained in the observations at different limb angles to simultaneously derive the oxygen 
abundance, A(O), and the scaling factor Sh that describes the cross-sections for inelastic collisions with neutral hydrogen relative the 
classical Drawin formula. Using the same codes and methods, we compare our 3D results with those obtained from the semi-empirical 
Holweger-Miiller model atmosphere as well as from different one-dimensional (ID) reference models. 

Results. With the CO 5 BOLD 3D solar model, the best fit of the center-to-limb variation of the triplet lines is obtained when the 
collisions by neutral hydrogen atoms are assumed to be efficient, i.e., when the scaling factor Sh is between 1.2 and 1.8, depending 
on the choice of the observed spectrum and the triplet component used in the analysis. The line profile fits achieved with standard ID 
model atmospheres (with fixed microturbulence, independent of disk position p) are clearly of inferior quality compared to the 3D 
case, and give the best match to the observations when ignoring collisions with neutral hydrogen (5h=0). The results derived with the 
Holweger-Miiller model are intermediate between 3D and standard ID. 

Conclusions. The analysis of various observations of the triplet lines with different methods yields oxygen abundance values (on 
a logarithmic scale where A(H)=12) that fall in the range 8.74 < A(O) < 8.78, and our best estimate of the 3D non-LTE solar 
oxygen abundance is A(O) = 8.76 ± 0.02. All ID non-LTE models give much lower oxygen abundances, by up to -0.15 dex. This is 
mainly a consequence of the assumption of a /^-independent microturbulence. An independent determination of the relevant collisional 
cross-sections is essential to substantially improve the accuracy of the oxygen abundance derived from the O i IR triplet. 

Key words. Sun: abundances, Sun: photosphere, Hydrodynamics, Radiative transfer, Line: profiles 


1. Introduction 

Oxygen, the most abundant chemical element of stellar nucle- 
osynthetic origin, is a widely used tracer for the evolution of 
various stellar populations. Oxygen is mainly produced in mas¬ 
sive stars that end their evolution as Type II supernovae. Since 
the main production of galactic iron by type la supernovae from 
less massive progenitor stars only contributes after some de¬ 
lay, the oxygen-to-iron abundance ratio, [O/Fefl can be used 
as an indicator of the chemical evolution of a particular com¬ 
posite stellar population. This abundance ratio, along with other 
a-element abundances, has been used to investigate Galactic 
star formation rates and to constrain models of the Galactic 
Chemical Evolution JMatteucci & Francoislil992i : lGratton et ali 


l2000i) . The solar oxygen abundance, AfOfl is a natural refer¬ 
ence for these investigations. 

Since oxygen is a volatile element, it condenses only partly 
in meteoritic matter, so that the relative amount of oxygen 
present in meteorites (e.g., with respect to a reference ele¬ 
ment like silicon) is not representative of the oxygen abun¬ 
dance ratio that was actually present in the original solar neb¬ 
ula. But it is thought that the latter is preserved in the so¬ 
lar photosphere, and hence can be inferred by spectroscopic 
analysis. Recent state-of-the-art solar spectroscopic abundance 
determinations resulted in a_rather low oxygen abundance of 
A(O) — 8.66 - 8.68 dex dPereira et al.1 l2009hl) . Previous esti¬ 
mates by Asnl und et al.l (120041) . who analyzed atomic and molec¬ 
ular oxygen lines in disk-integrated spectra, yielded a consistent 
value, A(O) = 8.66 + 0.05 dex. The work of lAllende Prieto et al.l 


1 [O/Fe] = log[N(0)/N(Fe)]* - log[N(0)/N(Fe)] 0 . 2 A(O) = log[N(0)/N(H)] + 12 
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( 2001 1 who analyzed the forbidden oxygen line at 630 nm, re- 
sulted in A(O) = 8.69 ± 0.05 dex, while the investigation of 
lAllende Prieto et aD (120041) . who performed a x 2 analysis of the 
3D non-LTE center-to-limb variation of the O IIR triplet, sug¬ 
gested a slightly higher value of A(O) = 8.72 dex. 


The rather low value of the solar oxygen abundance, A(O) Si 
8.70dex, as found in the above mentioned investigations, poses 
significant problems for theoretical solar structure models to 
meet the helioseismic constraints. These models depend on the 
solar chemical composition via radiative opacities, equation of 
state, and energy production rates. Most notably, changes in 
CNO abundances influence the structure of the radiative inte¬ 
rior, the location of the base of convective envelope, the he¬ 
lium abundance in the convection zone, and the sound speed 
profile. It was shown by Basu&Antia ( 20081) th a t the ‘ old’ so¬ 
lar chemical composition of Grev esse & S auvall d 19981) . where 
A(O) = 8.83+0.06, yields a much better agreement between the¬ 
ory and helioseismic observations than the new solar metallicity 
with a significantly reduced oxygen abundance of A(O) ^ 8.70. 


According to our own previous studies dCaffau et al.ll2008l 
hereafter Paper I), the solar oxygen abundance is somewhat 
higher, A(O) = 8.73 - 8.79 dex, with a weighed mean of 
A(O) = 8.76 ± 0.07 dex. This rather high oxyg en abundance 
was supported by the work of 1 Ayres et al.l d2013l) . who used a 
C05B0LD model atmosphere to analyze molecular CO lines (as¬ 
suming n(C)/n(Q)=0 .5) and derived A(O) = 8.78 + 0.02dex. 
IVillante et al.l (12014l) has shown that the chemical composition 
presented bv lcaffau et alJ d201 ill , which includes a ‘high’ A(O), 
is in much better agreement with the overall metal abundance 
inferred from helioseismology and solar neutrino data. 


Spectroscopic abundance analysis relies on stellar atmo¬ 
sphere models and spectrum synthesis, as a theoretical basis. 
The reliability of these abundance determinations depends on 
the physical realism of the atmosphere models and of the spec¬ 
tral line formation theory. On the observational side, high signal- 
to-noise spectra of sufficient spectral resolution are essential, as 
is the availability of clean spectral fines. The number of suit¬ 
able oxygen fines in the optical solar spectrum is very limited. 
The O i IR triplet lines at 777 nm are essentially free from any 
contamination, and hence could be expected to be a good abun¬ 
dance indicator. It is well known, however, that the triplet lines 
are prone to departures from local thermodynamic equilibrium 
(LTE) in solar-type stars (e.g.. lKiselmanl[l99ll) . The non-LTE 
effects lead to stronger synthetic spectral lines than in LTE, and 
hence a lower oxygen abundance is required to match the ob¬ 
served line strength. 


It has been repeatedly demonstrated that it is impossible to 
reproduce the center-to-l imb variation of the triplet lines unde r 
the assumpti on of LTE dKiselmanll 1 9931: IPereira et alJl20 09b). 
iPereira et alJ (l2009bl) have also shown that non-LTE synthetic 
fine profiles based on standard one-dimensional (ID) model at¬ 
mospheres (e.g., ATLAS, MARCS) cannot explain the observed 
properties of t he triplet, while the semi-empirical Holweg er- 
Miiller model dHolwegeil 1 1 9671 lHolweger~ & Muellei] Il974h is 
able to achieve a consistent match of the observations. Three- 
dimensional (3D) model atmospheres, however, give the best re¬ 
sults. 


In this paper, we expand on the work done in Paper I, fol¬ 
lowing up on the investigation of the solar Oi IR triplet at 
777 nm with an independent 3D, non-LTE analysis of several 
high-quality solar intensity spectra for different values of the di¬ 


rection cosine /0 (hereafter //-angle for short) and of a spectrum 
of the disk-integrated flux. In Section [3 we describe the spec¬ 
tra that we use for our analysis, and in Sect.[3] we outline the 
theoretical foundations of interpreting the spectra. The main re¬ 
sults are presented in Sect. [4] in which we describe the adopted 
methodology, perform different comparisons between observed 
and theoretical spectra, and derive the most likely value of the 
solar oxygen abundance. We discuss the results and their impli¬ 
cations in Sect. [3 before concluding the paper with a summary 
of our findings in Sect. [6] 


2. Observations 

We utilize intensity spectra of the O i IR triplet recorded at dif¬ 
ferent //-angles across the solar disk. As mentioned above, the 
analysis of the triplet lines requires a detailed modeling of their 
departures from LTE, which is rather sensitive to the strength of 
collisions with neutral hydrogen atoms. Since the relevant col- 
lisional cross-sections are essentially unknown, we derive both 
the oxygen abundance, A(O), and the scaling factor with re¬ 
spect to the classical (Drawin) cross-sections for collisional ex¬ 
citation and ionization of oxygen atoms by neutral hydrogen, 
.S'h, from a simultaneous fit of the line profiles or their equiva¬ 
lent widths at different //-angles. An analysis of the disk-center 
or disk-integrated spectrum alone would be insufficient to yield 
unique results. Since the impact of S'h depends systematically 
on the depth of line formation and therefore on //, the variation 
of the spectral lines with disk position adds another dimension 
to the analysis and helps to constrain the problem. Following a 
similar idea. ICaffau et all ( 2015 ) utilized the center-to-limb vari¬ 
ation of the oxygen and nickel blend at 630 nm to separate the 
contribution of the individual elements. 

We make use of the following three independent observa¬ 
tional data sets: 

1. The absolutely calibrated FTS spectra obtained at Kitt Peak 
in the 1980s, as described by INeckel & Labsl (1 19841) and 
iNeckelld 19991) . covering the wavelength range from 330 nm 
to 1250 nm for the disk center and full disk (Sun-as-a-star), 
respectively. In the following, we refer to this as “Neckel 
Intensity” and “Neckel Flux”. The spectral purity (A/l) 
ranges from 0.4pm at 330nm to 2pm at 1250nm (825 000 
> /1/A A > 625 000). 

2. A set of spectra observed by one of the authors (WCL) at 
Kitt Peak with the McMath-Pierce Solar Telescope and rapid 
scan double-pass spectrometer (see iLivingston et al.l 120071 
and references therein). Observations were taken in single¬ 
pass mode on 12 September, 2006; the local time was about 
1 Oh 15m to 10h45m. The resolution was set by the slit width: 
0.1 mm for the entrance and 0.2 mm for the exit slit. This 
yields a resolution of A/AA ~ 90000. The image scale is 
275 mm -1 . The slit has a length of 10mm and was set par¬ 
allel to the limb. Observations were made at six //-angles: 
// = 1.00, 0.87, 0.66, 0.48, 0.35, and 0.25, along the merid¬ 
ian N to S in telescope (geographic) coordinates. We have 
acquired two spectra for any // value and four spectra at disk 
center. The exposure time for each individual spectrum is 
100 s, thus essentially freezing the five-minute solar oscil¬ 
lations. No wavelength calibration source was observed, and 

3 // = cos 6, where 0 is the angle between the line of sight and the 
direction perpendicular to the solar surface; // = 1 at disk center, while 
// = 0 at the limb. 
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we used the available telluric lines to perform the wavelength 
calibration. The positions of the telluric lines were measured 
on the Neckel Intensity spectrum and a linear dispersion re¬ 
lation was assumed. 

After the wavelength calibration, we summed the spectra 
at equal /./ to improve the signal-to-noise ratio (S/N). The 
continuum of the summed spectra was defined by a cubic 
spline passing through a predefined set of line-free spectral 
windows. A comparison of the normalized disk-center 
spectrum with the Neckel Intensity atlas reveals a systematic 
difference in the equivalent widths, in the sense that the 
spectral lines are weaker in our own spectrum by about 
4%. We attribute this difference to instrumental stray light, 
and correct the normalized spectra for all //-angles as 
fcon-(d) = (1 + /) /(/!) - /, where / = 0.0405 is the stray 
light fraction. We estimate the signal-to-noise ratio of the 
stray light corrected spectra to be S/N ~ 1500 at disk center, 
and S/N » 700 at p — 0.25. In the following, we refer to 
these spectra as WCLC. 



3. The high-quality spatially and te mporally averaged inten¬ 
sity spectra used by Pere ira et al.1 (j2009b'FL recorded with 
the Swedish 1-m Solar Telescope (Scharmeret al. 20031) on 
Roque de Los Muchachos at La Palma over two weeks in 
May 2007. For these observations, the TRIPPEL spectro¬ 
graph was used with a slit width set to 25 pm, corresponding 
to 0.11" on the sky. The spectral resolution of this data set is 
/1/A/l « 200000. The data set consists of spectra at five disk 
positions, p = 1.00, 0.816, 0.608, 0.424, and 0.197. Since 
a large number of spectra at different instances of time and 
space were summed up, a high signal-to-noise ratio of S/N 
~ 1 200 - 1 500 was achieved at all //-angles. The spectra 
are corrected for stray light. See iPereira et al.l (i2009alb! ) for 
a detailed discussion regarding this observational data set. 

3. Stellar atmospheres and spectrum synthesis 

To compare observed and theoretical spectra, we performed de¬ 
tailed non-LTE spectral line profile calculations. In this section, 
we describe the theoretical basis of the calculations performed 
for the present analysis of the center-to-limb variation of the dif¬ 
ferent solar spectra. 


l0 9 T Ross 

Fig. 1. Top: temperature structure of 3D, (3D), ID LHD, and HM 
model atmosphere on the Rosseland optical depth scale. Bottom: 
horizontal temperature fluctuations, ATrms, of the 3D model 
at constant tr oss (black) and temperature difference HM - (3D) 
(dashed) and ID LHD- (3D) (blue) as a function of logrR 0ss . 


below, which controls the outgoing bolometric radiative flux (ef¬ 
fective temperature) of the model atmosphere. At the top bound¬ 
ary, the temperature of inflowing matter and the scale heights 
of pressure and optical depth beyond the upper boundary are 
controlled via free parameters. In the lateral directions, periodic 
boundary conditions are applied. Under these conditions, matter 
and radiation that exit the model nonvertically on one side, enter 
the box from the other side. 

The construction of the opacity tables that are needed for 
the radiat ion transport relie s on monochromatic MARCS opaci¬ 
ties (|Gustafsson et al. 20081), c omputed in LTE and adopti ng the 
solar chemical composition of iGrevesse & Sauvall (1 9981) . with 
the exception of the CNO abundances that were set to those 
of lAsplund et al.l d2005i) . To speed up the calculations, opac- 


3.1. Model atmospheres 

We give an overview of the properties of the 3D hydrodynamic 
and ID hydrostatic model atmospheres used in this paper. 


3.1.1.3D hydrodynamical C05B0LD model atmospheres 

The 3D hydrodynamical solar model atmosphere utilized in this 
work is very similar to that_used in Paper Q It was computed 
with the CO 5 BOLD code Frevtag et al.l ( 2012 ) on a Cartesian 
grid consisting of 140 x 140 grid points horizontally and 150 
grid points vertically, spanning 5.6 x 5.6 Mm 2 horizontally and 
2.25 Mm vertically (“box-in-a-star” setup). Open boundary con¬ 
ditions were applied vertically, allowing matter and radiation to 
flow in and out of the computational domain while conserving 
the total mass of matter inside the box. At the lower boundary, 
we prescribe the entropy of the matter entering the model from 


iti es are grouped in to 12 opacity bins accordin g to the ideas 
of Nordlund ( 1982b, an d furt her elaboration by Ludwi g et al.l 
(119941) and Vogler et al.l (120041) . In the calculation of the tables, 
the continuum scattering opacity was added to the true absorp¬ 
tion opacity. 

This setup was evolved in time until a thermally relaxed, sta¬ 
tistically steady state was reached. Subsequently, the run was 
continued for another two hours of stellar time, covering about 
24 convective turnover times. Out of this two-hour time se¬ 
quence, 20 representative model snapshots were chosen so that 
their mean effective temperature (radiative flux) and its tempo¬ 
ral RMS fluctuation closely agree with the corresponding val¬ 
ues of the fully resolved time sequence. In addition, we made 
sure that the horizontally averaged vertical velocity showed a 
similar depth dependence and temporal fluctuation in the sub- 
samp l e as in the whole sam ple. It was argued previously dCaffaul 
2009; iLudwig et al.ll2009h that such subsample of model atmo¬ 
spheres represents the whole evolution in terms of spectral line 


4 http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/5S8/148jroperties reasonably well. By using the selected snapshots, we 

5 The only difference is that we now employ an equation of state that can thus considerably reduce the computation cost of the subse- 


is consistent with the chemical abundance mix of the opacity tables. 


quent spectrum synthesis. For additional information about the 
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3D solar model used in this paper and comparison with other 
3D-hydrodynamical solar model atmospheres, see iBeeck et all 

d20l2) . 

We obtain the so-called (3D) model via horizontal and tem¬ 
poral averaging of temperature (T 4 ) and gas pressure (P) on sur¬ 
faces of constant optical depth (tr OS s)- The (3D) model is a ID 
atmosphere representing the mean properties of the 3D model 
(cf. Fig.Q}. Comparison of the line formation properties of the 
full 3D and the (3D) model gives an indication of the impor¬ 
tance of horizontal inhomogeneities in the process of line for¬ 
mation. In the present work, the performance of the (3D) model 
is compared with that of the standard ID models. 


package Linfor3D that was used to compute spectral line pro¬ 
files across the solar disk. In this section, we describe the basic 
assumptions and methods applied in the non-LTE spectrum syn¬ 
thesis calculations. 

3.2.1. The oxygen model atom 

We use an oxygen model atom that consists of 22 levels of O i 
(plus continuum), each having an associated bound-free transi¬ 
tion. In addition, we take 64 bound-bound transitions into ac¬ 
count, some of which are purely collisional. Fig. [2] shows the 
Grotrian diagram of our model atom. 


3.1.2. 1D LHD model atmospheres 

In the present analysis, we have also used classical ID model at¬ 
mospheres computed with the LHD code (see lCaffau et alJ2008lh 
We computed ID LHD and 3D C05B0LD model atmospheres us¬ 
ing identical stellar parameters, chemical composition, equation 
of state, opacity binning schemes, and radiation transport meth¬ 
ods, and thus these models are suitable for the analysis of 3D-1D 
effects through a differential comparison. Since LHD model at¬ 
mospheres ar e static, convection is treated via the mixing length 
theory (MLT. lB5hm-Vitensell 1958b iMihalasll 1 9781) . 

Since the oxygen triplet lines have a high excitation poten¬ 
tial and, hence, form in the deep photosphere, it may be expected 
that their strength is sensitive to the choice of mixing length pa¬ 
rameter, q;mlt- To probe the possible influence of o-yilt, we have 
used a grid of LHD model atmospheres with ctmlt = 0.5,1.0,1.5. 
After performing some tests, however, we found negligible in¬ 
fluence of otmct on the triplet lines (see Appendix[D]for details). 
Hence, all results shown in the following are based in the LHD 
model atmosphere with o'mlt = 1 -0. 


3.1.3. Semi-empirical Holweger-Muller solar model 
atmosphere 

W e have also used the semi-em pirical model atmosphere derived 
bv lHolweger & Muellerldl 9741) (in the following HM model, see 
also Holweger 1967). The T(t) relation of this model was em¬ 
pirically fine-tuned to match observations of the solar continuum 
radiation and a variety of spectral lines. 

The temperature structure of the HM model is significantly 
different from that of the theoretical ID LHD model atmo¬ 
spheres, as illustrated in Fig.Q] In the optical depth range -5.0 < 
logTR OS < -0.5, the HM model is slightly warmer than the 
(3D) model by 50 to 100 K, while the ID LHD model almost 
perfectly matches the (3D) model. In the deep photosphere 
(-0.5 < log tr os < +0.5), where significant parts of the triplet 
lines originate, the HM model closely follows the temperature 
of the (3D) model, while the ID LHD shows a somewhat steeper 
temperature gradient in these layers. This results in stronger OI 
triplet lines computed from the LHD model compared to those of 
the HM atmosphere. 


3.2. 3D/1D non-LTE spectrum synthesis 

Theoretical spectral lines were synthesized in two separate steps. 
Firstly, the departure coefficient^ of the atomic levels involved 
in the O i triplet were computed with the NLTE3D code. Then the 
departure coefficients were passed to the spectral line synthesis 



Fig. 2. Oxygen model atom used in this work. Radiatively al¬ 
lowed transitions are marked with green solid lines. Radiatively 
forbidden transitions, which were treated via electron-impact ex¬ 
citation only, are marked by blue dashed lines. The forbidden 
[O i] line and the O i IR triplet are marked by red dash-dotted 
lines. We took collisional and radiative bound-free transitions 
into account for each level. 


Level energies and statistical weights of the model atom 
were taken from the NIST Atomic Spectra Database Version fQ| 
(see Tab. lA. Ik Cross sections for radiative, bound-free transi¬ 
tions were taken from the TOPBASE servioj] of the Opacity 
Project (ICunto et al.l 119931) . Lacking reliable collisional cross 
sections derived from laboratory measurements or quantum- 
mechanical cal culations, w e resort to the clas sical prescrip¬ 
tions of Allen (fAl lenl 1 197 3L p. 41) and Drawin (lDrawinlll969l : 
ISteenbock & Holwegedll984i lLamberdll993h to obtain the rate 
coefficients for bound-free transitions due to ionizing collisions 
with electrons and neutral hydrogen, respectively. As the Drawin 
approximation is highly uncertain (see Sect. l5.3l bclow). we scale 


6 The departure coefficient of atomic level i, b n is defined as the ratio 6 7 http://www.nist. gov/pml/data/asd. cfm/ 

of non-LTE to LTE population number density, n f /n* 8 http ://cdsweb .u-strasbg. fr/topbase/topbase .html 
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the rates obtained from the Drawin formula by the adjustable pa¬ 
rameter .S'h , which is deduced simultaneously with the oxygen 
abundance (see Sect.[4}. Finally, we include the charge transfer 
reaction O l) + H + ^ C)_ + H° for th e grou nd level of oxygen, 
as described bv lArnaud & Rothenflud d 1985b . Because the ion¬ 
ization potentials of O i and HI are similar, this reaction is nearly 
resonant and very efficient. It ensures that the ground level of O i 
is in equilibrium with O n. 

Out of the 64 bound-bound transitions, 54 are radiatively 
and collisionally permitted, while the remaining ten are radia¬ 
tively forbidden. For the latter transitions, we only take colli¬ 
sions with electrons into account (see below). Einstein coeffi¬ 
cients and wavelengths for the permitted transitions were again 
taken from the NIST Atomic Spectra Database (see Tab. lA.2l >. 

Rate coefficients for collisional excitation by neutral hydro¬ 
gen were computed with classical Drawin formula and scaled by 
the same factor S'h as applied to the ionizing collisions by neu¬ 
tral hydrogen. Since the collisional cross section according to 
the Drawin formalism is proportional to the oscillator strength 
gf of the transition, the levels of the radiatively forbidden tran¬ 
sitions are also not coupled by hydrogen collisional excitation 
in our model atom. Alternative assumptions and their impact on 
the derived oxygen abundance are briefly discussed in Sect. 15.3.11 
below. 

Rate coefficients for the collisional excitation by e~ w ere 
taken from the detailed R-matrix computations of iBarkleml 
d2007l) . For the radiatively allowed transitions that were not 
included in the latter work, we applied the classical van 
Regemor ter’s formula (lvan~R egemorten ll962l) . The rate coeffi¬ 
cients of lBarkleml (l2007h . unlike those computed by the oscilla¬ 
tor strength-dependent classical van Regemorter’s formula, are 
nonzero for radiatively forbidden transitions (blue dashed lines 
in Fig.O, and hence we were able to account for the inter-system 
collisional transitions. The latter transitions are thought to play 
an important role in defining the statistical equilibrium of oxy¬ 
gen (especially at low metallicities, see lFabbian et al.li2009h and 
must be taken int o ac count in realistic NLTE simulations. While 
the work of iBarkleml d2007l) provides data for a total of 171 
bound-bound transitions, we have only used 48 of them, based 
on the magnitude of their rate coefficient. We have performed 
tests that evaluate the influence of the remaining purely colli¬ 
sional transitions and found that it is negligible in the case of the 
Sun. 

It should be noted that fine-structure splitting is ignored in 
our model atom, and all multiplets are treated as singlets. In 
particular, the upper level of the Oi IR triplet, level 3p 3 P, is 
treated as a single superlevel with a single departure coefficient. 
However, since the three triplet components are well separated 
in wavelength, the line opacity computed from the combined gf- 
value of the superlevel transition would be too large, resulting in 
a wrong mean radiation field and, in turn, in wrong radiative ex¬ 
citation rates. To avoid this problem, we replace - for the radia¬ 
tive line transfer calculations only - the summed-up g/-value of 
the combined multiplet transition with the g/-value of the inter¬ 
mediate triplet component, thus diminishing the line opacity by 
a factor 3 (see Appendix iBl for additional information). Similar 
opacity scaling is applied for all other lines that involve multi¬ 
plet transitions. This approach significantly reduces the compu¬ 
tational cost of solving the statistical equilibrium equations. 

3.2.2. Computation of departure coefficients with NLTE3D 

As a first step towards computing 3D/1D non-LTE line profiles 
of the O i IR triplet, we computed departure coefficients bj, using 


the NLTE3D code with the 3D/1D solar model atmospheres and 
oxygen model atom described in the previous sections. 

Firstly, NLTE 3D computes the line-blanketed radiation field 
that is required for the calculation of the photoionization transi¬ 
tion probabilities. For th is task, we utilized the li ne opacity dis¬ 
tribution functions from ICastelli & Kuruczi d2004 ) and the con¬ 
tinuum opacities from the IONOPA package. The radiation field 
is derived with a long-characteristics method where the trans¬ 
fer equation is solved with a modified Feautrier scheme. To save 
computing time, continuum scattering is treated as true absorp¬ 
tion. For solar conditions, the validity of the latter approxima¬ 
tion has been verified via comparison with test calculations for a 
single 3D snapshot in which coherent isotropic continuum scat¬ 
tering is treated consistently. To further minimize the compu¬ 
tational cost, we have reduced the spatial resolution of the 3D 
model by using only every third point in both horizontal direc¬ 
tions, thus saving about 90% of the computing time. Tests have 
shown that this loss of spatial resolution does not influence the 
final results in any significant way (see Appendix 0. 

Based on this fixed 3D radiation field, we compute the 
photoionization and photorecombination transition probabilities. 
They are assumed to be independent of possible departures from 
LTE in the level population of O i. Then, NLTE3D computes the 
angle-averaged continuum intensities, / cont , at the central wave¬ 
length of all of the line transitions. If a particular spectral line 
is weak, the profile averaged mean line intensity, J is fixed to 
■/cont, and the radiative bound-bound transition probabilities of 
these lines are only computed once. Finally, we compute the 
collisional transition probabilities, which depend only on the lo¬ 
cal temperature. All of those transition probabilities are fixed, as 
they do not depend on the departure coefficients. In this situa¬ 
tion, the solution of the statistical equilibrium equations would 
be straightforward, and the departure coefficients of all levels 
could be obtained immediately, i.e., without any iterative proce¬ 
dure. 

However, as soon as the line transitions are no longer opti¬ 
cally thin, the mean radiative line intensities, /, and the depar¬ 
ture coefficients, bj, are inter-dependent quantities that exhibit a 
highly non-linear coupling. In this case, we must compute / as 
the weighted average of J v over the local absorption profile, tak¬ 
ing into account that line opacity and source function depend on 
the departure coefficients of the energy levels of the transition. 
In the present version of NLTE 3D, the local line profile is as¬ 
sumed to be Gaussian, representing the thermal line broadening 
plus the effect of microturbulence. At this point, the 3D hydro- 
dynamical velocity field is not used, so Doppler shifts along the 
line of sight are ignored. Instead, the velocity field is represented 
by a depth-independent microturbulence of 1.0 km/s. Radiative 
damping, as well as van der Waals and Stark broadening, are 
ignored at this stage. We verified that this simplification has a 
negligible influence on the resulting departure coefficients. For 
the spectrum synthesis, however, the line broadening is treated 
in full detail and correctly includes the effects of the 3D hydro- 
dynamical velocity field (see below). 

Since it would be very computationally expensive to 
use complete linearization methods (lAuerl[l973t iKubat 2010, 
Sect. 8) with 3D model atmospheres, we choose to iterate be¬ 
tween mean radiative line intensities and depa r ture c oefficients. 
However, as discussed by Prakanavicius et al. l ( 2013 ). the solu¬ 
tion of the oxygen non-LTE problem with NLTE3D by ordinary 
A-iteration exhibits a poor convergence, which could potentially 
lead to faulty results. Hen ce, we implemented the acce lerated A- 
iteration (ALI) scheme of iRvbicki & Hummed (1199ll Sect. 2.2) 
to achieve faster convergence. 
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Iterations are started with LTE line opacities and source 
functions, from which the mean line intensities, /, radiative tran¬ 
sition probabilities, and equivalent widths (EWs) of the emer¬ 
gent flux line profiles are calculated. Then, all previously com¬ 
puted transition probabilities enter the equations of statistical 
equilibrium, which are solved for the unknown departure coeffi¬ 
cients. Using the computed departure coefficients, the non-LTE 
line opacities and source functions are computed and the radia¬ 
tion transfer is solved for all relevant lines using these quantities. 
Then the J, flux EWs, and approximate A-operators (A*) are re¬ 
computed with non-LTE line opacities and source functions. The 
updated mean radiative intensities enter the equations of statisti¬ 
cal equilibrium and the cycle begins again. 

The accelerated A-iterations are continued until the relative 
changes (from one iteration to the next) in the flux EWs of all 
sufficiently strong lines have converged to less than 10 5 . We 
checked that increasing the accuracy of the convergence crite¬ 
rion does not alter the resulting departure coefficients signifi¬ 
cantly, which might be ascribed to the good convergence prop¬ 
erties of the accelerated A-iteration scheme. 

In preparation for the following investigation, we computed 
departure coefficients for a grid of Sh and A(O) values. The .S'h 
values range from 0 to 8/3 in steps of 1/3, and A(O) varies from 
8.65 to 8.83 dex in steps of 0.03 dex. The spectrum synthesis cal¬ 
culations described in Sect. l3.2.3l are based on this rather exten¬ 
sive and dense grid of departure coefficients. 

Since the NLTE3D code can only work with 3D model at¬ 
mospheres defined on a geometrical grid, the ID LHD and HM 
model atmospheres were converted to pseudo-3D models by 
replicating the particular ID model atmosphere three times hori¬ 
zontally in each direction so that the resulting model would have 
3x3 identical columns. Using the same equation of state and 
Rosseland opacity tables as in the 3D models, the vertical grid 
spacing of the pseudo-3D models is adjusted to ensure that the 
resulting model atmosphere is sampled at equidistant log tr oss 
levels, and preserves the T(t) relation of the original ID model. 


3.2.3. 3D/1D NLTE spectrum synthesis with Linfor3D 

The second step in computing the synthetic line profiles of the 
oxygen infrared triplet is the spectrum synthesis itself. Here we 
adopt the same atomic line parameters as specified in Table 1 of 
Paper I, and perform detailed 3D non-LTE line formation calcu¬ 
lations with the Linfor3D0 package, which can handle both 3D 
and ID model atmospheres. As in the case of the NLTE3D calcu¬ 
lations of the departure coefficients, we reduced the horizontal 
resolution of the 3D model atmosphere by choosing only every 
third grid point in x and //-direction. Again, we conducted exper¬ 
iments confirming that this reduction does not influence the final 
results (see also Appendix O. 

Linfor3D computes the LTE populations of the upper 
and lower levels involved in the triplet transitions via Saha- 
Boltzmann statistics using the IONOPA package^- From this, 
Linfor3D evaluates the LTE continuum and line opacities, /r/ljj 1 . 
and /rj'j' L K . Then it uses the previously derived departure coeffi- 


9 http://www.aip.de/Members/msteffen/linfor3d 

10 Since this package is also utilized in NLTE3D, the calculations of de¬ 
parture coefficients and spectrum synthesis share consistent LTE num¬ 
ber densities and continuum opacities. 


cients to compute the NLTE line opacity and line source 
function S l/j/' rK (assuming complete redistribution), according to 

4lte , exp i%)-{bjb0 

—r— -b\ - i - (1) 

«LTE exp(&)-l 

and 

Snlte = exp(&)-l 

By (bi/b u ) exp(^) - 1 ’ 

where B v is the Planck function, and b\ and b u are the departure 
coefficients of the lower and up per le vel of the line transition, 
respectively (see, e.g., lMihalaslll978l . Chapter 4-1). Note that 
both LTE and NLTE line opacities are corrected for stimulated 
emission. We recall that we assume all transitions of the oxy¬ 
gen 777 nm triplet to share identical departure coefficients b\ and 
b u . The local line profiles are now represented by Voigt profiles, 
taking into account thermal broadening combined with van der 
Waals and natural line broadening]; differential Doppler shifts 
along the line of sight due to the 3D hydrodynamical velocity 
field are considered as well. 

Using these ingredients, the radiative transfer equation is in¬ 
tegrated on long characteristics to obtain the emergent radiative 
intensities as a function of wavelength. The intensity spectra are 
calculated for two sets of inclination angles //, which represent 
the WCLC and iPereira et al.l (12009b ) observations (see Sect.|2] 
or Table 0, in each case for four azimuthal directions (0, 7r/2, n, 
3/2 n). The final intensity spectrum for each p is obtained by av¬ 
eraging over all surface elements, and to improve the statistics, 
over all azimuthal directions. For comparison with the solar flux 
spectrum, we computed the integral over p of the p-weighted in¬ 
tensity spectra using cubic splines to derive the disk-integrated 
flux profile. 

We finally compute a grid of synthetic line profiles using the 
corresponding grid of departure coefficients described in the pre¬ 
vious section. Each component of the oxygen triplet is synthe¬ 
sized for .S’h- values in the range from 0 to 8/3 (with steps of 1/3), 
and A(0)-values in the range between 8.65 and 8.83 (with steps 
of 0.03 dex). 


4. Results 

Computed synthetic intensity spectra are compared with obser¬ 
vations across the solar disk to derive A(O) and .S’h simultane¬ 
ously. For this task, we adopt two different methodologies in 
which we compare theoretical and observed line profiles and line 
equivalent widths, respectively. In this section, we describe the 
fitting procedures and their individual results. 

4. 1. Line profile fitting 

Here we compare observed and theoretical intensity line pro¬ 
files and search for the best-fitting A(0)-Sh combination. This 
task utilizes the grid of theoretical NLTE line profiles described 
in Sect. 13.2.31 ) and searches for the A(0)-.S)| combination that 
provides the best fit of the intensity profiles of a single triplet 
component simultaneously at all //-angles. 


11 Stark broadening is ignored here as in Paper I. 
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Fig. 3. Best-fitting non-LTE line profiles computed with the 3D model atmosphere for the WCLC (top) and Pereira (bottom) spectra. 
The top panels show the observed spectrum (black lines) and their theoretical best-fitting counterparts (red dots) for the different 
p-angles, p increasing from top (limb) to bottom (disk center). For clarity, a vertical offset of 0.1 was applied between consecutive 
p-angles. The legend shows the best-fitting A(O) value for each of the three triplet lines. The bottom panels show the difference of 
the normalized intensites 7 obs - 7 calc . An offset of 0.01 was applied between the different //-angles (ft increasing downwards). 


4.1.1. Fitting procedure 


We employ the Levenberg-Markwardt least-squares fitting al¬ 
gorithm (iMarkwardd 120091) . as implemented in 1DL procedure 
MPFIT, to find the minimum y 2 in the four-dimensional pa¬ 
rameter space defined by A(O), Sh, i/b, and Ad. Here iy is the 
FWHM (in velocity space) of a Gaussian kernel, which is used 
to apply some extra line broadening to the synthetic spectra 
to account for instrumental broadening and, in the case of ID 
model atmospheres, macroturbulence; A/l is a global wavelength 
shift applied to the synthetic spectrum to compensate for poten¬ 
tial imperfections of the wavelength calibration, uncertainties in 
the laboratory wavelengths, and systematic component- and //- 
dependent line shifts. Test runs with an additional free fitting 
parameter for the continuum placement resulted in unsatisfac¬ 


tory fits. We therefore rely on the correct normalization of the 
observed spectra and consider the continuum level as fixed. 

In a first sweep, we find the minimum y 2 for each triplet 
component k and each /./-angle individually. Here y 2 is defined 
as 



where the sum is taken over all n considered wavelength points 
of the line profile, 7 obs and 7 ( calc is the intensity of observed and 
synthetic normalized spectrum, respectively, and cr, is the pixel- 
dependent observational uncertainty of the spectrum to be fit¬ 
ted. For simplicity, we have adopted a fixed continuum signal- 
to-noise ratio of S /N — 500 for both the WCLC and the Pereira 
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Table 1. Results of non-LTE line profile fitting with various model atmospheres: A(O), .S'n-values, and reduced^ 2 of the best fit 
for two different observed spectra (WCLC and Pereira 2009); the difference between A(O) derived from the two different spectra is 
given in Col. (4). For the individual triplet components, we give the formal fitting errors of A(O) and ,S'n (only if .S'h> 0), while the 
errors of the mean values correspond to the standard deviation of the individual lines. 



A(O) 

AA(O) 

Sh 

y 2 

Data set: 

WCLC 

Pereira 


WCLC 

Pereira 

WCLC 

Pereira 

d, nm 

3D model atmosphere 

777.2 

8.748 ± 0.002 

8.759 ± 0.004 

-0.011 

1.44 ±0.03 

1.64 ±0.07 

1.654 

1.794 

777.4 

8.735 ± 0.002 

8.769 ± 0.004 

-0.033 

1.22 ±0.03 

1.67 ±0.07 

1.073 

0.853 

777.5 

8.756 ± 0.002 

8.771 ±0.003 

-0.015 

1.58 ±0.04 

1.85 ±0.07 

1.556 

1.157 

Mean 

8.747 ±0.010 

8.766 ± 0.006 

-0.019 

1.41 ±0.19 

1.72 ±0.12 

- 

- 

d, nm 


(3D) model atmosphere, £ mic = 

l.Okm/s 



777.2 

8.591 

8.641 ±0.011 

-0.050 

0.00 

0.22 ± 0.05 

6.163 

7.010 

777.4 

8.603 

8.609 

-0.006 

0.00 

0.00 

3.985 

3.522 

777.5 

8.614 

8.623 ±0.010 

-0.009 

0.00 

0.03 ± 0.05 

2.718 

2.358 

Mean 

8.603 ±0.011 

8.624 ±0.016 

-0.021 

- 

0.08 ±0.12 

- 

- 

d, nm 


HM model atmosphere, ^nu C = 

0.8km/s 



777.2 

8.671 ±0.004 

8.673 ± 0.008 

-0.002 

0.44 ± 0.03 

0.46 ± 0.06 

4.409 

4.387 

777.4 

8.680 ± 0.004 

8.713 ±0.007 

-0.033 

0.43 ± 0.03 

0.66 ± 0.06 

3.042 

2.322 

777.5 

8.747 ± 0.002 

8.759 ± 0.005 

-0.012 

1.02 ±0.03 

1.18 ±0.07 

2.497 

2.006 

Mean 

8.699 ±0.041 

8.715 ±0.043 

-0.016 

0.63 ± 0.34 

0.76 ±0.37 

- 

- 

d, nm 


HM model atmosphere, = 

1.2km/s 



777.2 

8.635 ± 0.005 

8.638 ±0.010 

-0.003 

0.38 ±0.03 

0.40 ± 0.06 

4.935 

4.998 

777.4 

8.644 ± 0.005 

8.671 ±0.007 

-0.027 

0.34 ± 0.03 

0.49 ± 0.06 

3.349 

2.652 

777.5 

8.715 ±0.002 

8.727 ± 0.005 

-0.012 

0.86 ± 0.03 

1.00 ±0.07 

2.674 

2.189 

Mean 

8.664 ± 0.044 

8.679 ± 0.045 

-0.015 

0.53 ± 0.29 

0.63 ± 0.32 

- 

- 

d, nm 

LHD model atmosphere, a MLT = 1.0, £ mic = 0.8km/s 

777.2 

8.608 

8.607 

0.001 

0.0 

0.0 

12.477 

14.021 

777.4 

8.616 

8.620 

-0.004 

0.0 

0.0 

8.661 

8.090 

777.5 

8.623 

8.625 

-0.002 

0.0 

0.0 

4.973 

4.689 

Mean 

8.616 ±0.008 

8.617 ±0.010 

-0.001 

- 

- 

- 

- 

d, nm 

LHD model atmosphere, qtmlt = 1,0, £ m ic = 1.2km/s 

777.2 

8.581 

8.580 

0.001 

0.0 

0.0 

14.064 

15.881 

777.4 

8.593 

8.597 

-0.004 

0.0 

0.0 

9.737 

9.264 

777.5 

8.604 

8.606 

-0.002 

0.0 

0.0 

5.548 

5.310 

Mean 

8.593 ±0.011 

8.594 ±0.013 

-0.001 

- 

- 

- 

- 


spectra, and computed cr,- simply as cr = (S/N) 1 , independent 
of wavelength point i and inclination angle //. 

An obvious blend in the blue wing of the strongest triplet 
component at A 111 .2 nm is excluded from the relevant wave¬ 
length range by masking a window between -374 and -214 mA 
from the line center, setting 1/cr,- = 0 for these pixels. 

The result of this x 2 minimization procedure is a set of best- 
fit parameters A(O), .S’h, (- ! b, A/l for each (k.fj). These interme¬ 
diate results are only used for two purposes: (i) for fixing the 
relative wavelength shift of the line profiles as a function of // 
for each triplet component k\ and (ii) to ‘measure’ the equiva¬ 
lent width of each individual observed line profile via numerical 
integration of the best-fitting synthetic line profile over a fixed 
wavelength interval of +2 A. The measured equivalent widths 
are used later for the equivalent width fitting method as described 
in Sect. 14. 2. Tl 

We note that the oxygen abundance obtained from fitting the 
individual line profiles is relatively high, 8.78 < A(O) < 8.87, 
where the required oxygen abundance tends to increase towards 
the limb. At the same time, the corresponding .S'h- values are 
large, too (Sh A 3.0). However, we do not consider these results 
to be relevant because the line profile at a given single // can be 
reproduced almost equally well with different combinations of 
the fitting parameters A(O), S'h, and t/b, which are tightly cor¬ 
related. Small imperfections in the theoretical and / or observed 


line profiles can thus easily lead to ill-defined solutions. This de¬ 
generacy is only lifted when fitting the profiles at all //-angles 
simultaneously. 

In the second sweep, we fit the line profiles at all //-angles 
(jUj, j = 1... in) simultaneously, but separately for each of the 
three triplet components. In this case, x 1 is defined as 


A" 


m n 


EE 

;=i i=i 



(4) 


where m is the number of considered //-angles (jj Z 0.2). Again, 
cTj j is taken to be constant, <x,j = cr = (S /A0~', as explained 
above, and the blend in the blue wing of the strongest triplet 
line at A 111 .2 nn is masked as before. Treating the triplet com¬ 
ponents separately has the advantage of allowing a consistency 
check and providing an independent error estimate of the derived 
oxygen abundance. 

As in the first sweep, we consider A(O), Sh, and Vb as free fit¬ 
ting parameters, but instead of a single global Ad, we introduce 
two global line shift parameters, Adi and Ad 2 , to allow for the 
most general case of a /./-dependent line shift correction, com¬ 
puted as Ad(//) = Ado(/i) + Adi + Ad 2 (1 - //), where AAo(/u) is 
the individual wavelength shift determined in the first sweep for 
the respective //-angle. It turns out, however, that the results of 
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the global fitting are insensitive to the detailed treatment of the 
line shifts. 


4.1.2. Fitting results for the WCLC and Pereira spectra 

Fig. [3] shows the results of the simultaneous line profile fitting 
with the 3D-NLTE synthetic line profiles of the WCLC and 
Pereira spectrum, respectively. In general, the agreement be¬ 
tween observation and model prediction is very good; the dif¬ 
ferences never exceed 1% of the continuum intensity and usu¬ 
ally remain at a much lower level. Both the line wings and line 
core are represented very well at all inclination angles. The very 
satisfactory match we obtain between theory and observation 
suggests that both the velocity field and the thermal structure of 
our model atmosphere are realistic in the range of optical depths 
sampled by the O i IR triplet lines for inclination angles p k. 0.2. 
In addition, the accurate reproduction of the observed line pro¬ 
files suggests that our NLTE line formation technique is ade¬ 
quate, despite the possible shortcomings regarding the modeling 
of hydrogen collisions. 

The corresponding fits obtained with the ID hydrostatic 
models are collected in Appendix[E] In general, the line profiles 
computed with the ID hydrostatic model atmospheres result in a 
less satisfactory agreement with the observations. For example, 
given a single A(O) and Sh value, the (3D) and ID LHD mod¬ 
els fail to simultaneously account for the line core and wings 
at // = 1, and predict a much too weak line core close to the 
limb (see Figs. IE.H and IE. 3I >. The situation with the HM model 
(Fig . IE.2b is better, but the line profile fits are still clearly infe¬ 
rior to those computed with the 3D model atmosphere. It is worth 
mentioning that, for all ID model atmospheres, changing the (p- 
independent) microturbulence parameter only changes the value 
of the best-fitting A(O) and .S'n, but does not improve the quality 
of the fits. 

In Table [Q we summarize the main results derived from 
fitting the line profiles simultaneously for all //-angles, sepa¬ 
rately for each of the three triplet components, and for both 
the WCLC and the Pereira observed spectra. Columns (2) and 
(3) list the best-fitting oxygen abundance A(O) and their formal 
fitting errors for WCLC and Pereira, respectively, and Col. (4) 
shows the difference between A(O) derived from the two differ¬ 
ent sets of spectra. Columns (5) and (6) give the corresponding 
results for .S'h, while the last two columns show the reduced % 2 
of the best fit for the two sets of spectra, respectively, defined as 
Tred = X 1 K n x m ~ P)i where p = 5 is the number of free fitting 
parameters. The fact that is of the order of unity means that 
the fit provided by the synthetic spectra is good (at least for the 
3D model), given the assumed uncertainty of the observed spec¬ 
tral intensity of cr, = 1/500. However, only relative -values 
are considered meaningful forjudging the quality of a fit. 

The formal fitting errors of A(O) and .S'n given in Tableware 
computed as the square root of the corresponding diagonal ele¬ 
ments of the covariance matrix of the fitting parameters, which 
is provided by the MPFIT procedure. The fact that these fitting 
errors are relatively small shows that the solution found by the 
fitting procedure is well defined. However, these formal errors 
do not account for possible systematic errors (e.g., due to imper¬ 
fect continuum placement, unknown weak blends, insufficient 
time averaging, etc.), and hence the true errors are likely to be 
larger. In the case of the (3D) and LHD models, the formal errors 
are meaningless because the best fit would require Sh< 0 (see 
below), and the constraint ,S'n> 0 essentially enforces 5 h=0. 


As shown in Table [T| the oxygen abundances derived from 
the two sets of observed spectra generally agree within 0.03 dex 
or better, which is of the same order of magnitude as the line- 
to-line variation of A(O). The WCLC spectra produce systemat¬ 
ically lower abundances. We attribute this to the systematically 
stronger lines at low p compared to the iPereira et al.l d2009bi) 
data (cf. Fig.[4]below), which implies a reduced gradient W(p), 
and, in turn, a lower .S'n and A(O), as explained in more detail in 
Sect. 15.21 In the case of the 3D model atmosphere, we consider 
the agreement in the derived A(O) between the different triplet 
lines and between the two sets of spectra as satisfactory. 

In case of ID LHD model atmospheres, the agreement be¬ 
tween the best-fit parameters of the two spectra is nearly perfect, 
but this is again related to the fact that .S'n is forced to 0. The 
resulting fits are not good (high ^ ,), particularly at the lowest 
//-angle, as mentioned above. The derived A(O) are significantly 
lower (by about -0.15 dex) than for the 3D case. Similar con¬ 
clusions can be drawn for the (3D) model. The non-LTE oxygen 
abundance obtained with the ID semi-empirical HM model at¬ 
mosphere lies somewhere in between the 3D and the ID LHD 
results. While S'h is positive, the line profiles are poorly fitted, 
both at disk center and near the limb (see Fig. E3, and the line- 
to-line variation of A(O) is rather high. 

As will be discussed below in Sect. 15.21 all ID results are 
systematically biased towards low Sh and low A(O) values. In 
the following, we shall therefore ignore the ID results for the 
determination of the solar oxygen abundance. 

To estimate the solar oxygen abundance from the profile fit¬ 
ting method, we compute, separately for each of the two data 
sets, the average abundance, (A(O)), as the unweighted mean 
over the best-fitting A(O) values of the three triplet lines. These 
mean oxygen abundances are listed in Table Q] in the rows 
marked ‘Mean’, together with their errors, which in this case are 
defined by the standard deviation of the individual data points 
from the three triplet components. 

From the 3D-NLTE synthetic spectra we find: (A(O)) = 
8.747 ± 0.010dex (WCLC spectra) and (A(O)) = 8.766 + 
0.006 dex (Pereira’s spectra). We note that the 1 cr error bars of 
the two results fail (by a narrow margin) to overlap, indicative of 
unaccounted systematic errors. We finally derive a two-spectra 
average of «A(0)» = 8.757 ± 0.01 dex («S H » = 1.6 ± 0.2). 

4.1.3. Fitting of disk-center intensity and flux line profiles 

As an independent consistency check, we also applied our line 
profile fitting procedure to the lNeckel & Labsl (1 1984b FTS spec¬ 
tra. In this case, we simultaneously fit the solar disk-center 
(Neckel Intensity) and disk-integrated (Neckel Flux) Oi IR 
triplet line profiles. The blue wing of the 777.2 nm line is masked 
as before. The synthetic flux spectrum is computed by integra¬ 
tion over p of the //-weighted intensity spectra available at the 
five inclinations defined by the Pereira spectra. Subsequently, 
the flux spectrum is broadened with the solar synodic projected 
rotational velocity v sin i = 1.8km/s. The fitting parameters are 
again A(O), .S'h, instrumental broadening //b, and the individual 
line shifts Adi and AU for intensity and flux profile, respec¬ 
tively. 

Using the 3D non-LTE synthetic line profiles, we obtain an 
average oxygen abundance of A(O) = 8.80 ± 0.008, where this 
error refers to the formal fitting error; the line-to-line scatter is 
much smaller, F 0.001 dex. The best-fitting mean .S'n value is 
2.7 ± 0.25. As before, the quality of the fits is very satisfactory, 
with = 0.81, 0.84, and 1.13 for the three triplet components 
at A 777.2, 111 A, and 777.5 nm, respectively. 
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Both A(O) and .S’n are significantly higher than inferred from 
the center-to-limb variation of the line profiles in the WCLC 
and Pereira spectra. We have to conclude that neither the formal 
fitting errors nor the line-to-line variation of A(O) represent a 
valid estimate of the true uncertainty of the derived oxygen abun¬ 
dance. Rather, systematic errors must be the dominating factor. 
We shall demonstrate in the next Section that the result of simul¬ 
taneously fitting the intensity and flux profiles is very sensitive to 
small errors in measuring or modeling the relative line strengths 
in the disk-center and integrated disk spectra. We therefore con¬ 
sider the oxygen abundance derived from fitting the spectra at 
several //-angles simultaneously as more reliable. 


4.2. Equivalent width fitting 

Another possible methodology for estimating the solar oxygen 
abundance is based on the comparison of the observed equiva¬ 
lent widths (EWs) and their center-to-limb variation with their 
theoretical counterparts. Matching the EW of the lines is admit¬ 
tedly less precise than line profile fitting: it is possible that syn¬ 
thetic and observed EWs match perfectly, while the line profiles 
exhibit significantly different shapes. However, EW-fitting has 
the advantage of being simpler: there are only two free fitting 
parameters, A(O) and .S'h. The additional fitting parameters that 
were necessary in the case of line profile fitting, related to extra 
line broadening and line shifts, have no influence on the equiva¬ 
lent width and thus are irrelevant for the EW fitting procedure. 


Table 2. EWs of the Oi IR triplet as a function of g measured 
from the different spectra with our line profile fitting method 
(this work); for the Pereira spectra, we also give the original 
measurements taken from IPereira et aTJ (f2009bl) . 



Equivalent width 

(pm) 


Oi7772 

Oi7774 

Oi7775 

p 

Spectrum: Neckel & Labs (1984). EW: this work 

1.00 

8.73 ± 0.050 

7.60 ± 0.044 

6.01 ±0.035 

flux 

7.54 ± 0.043 

6.52 ±0.038 

5.07 ± 0.029 

P 

Spectrum: WCLC, EW: this work 

1.00 

8.73 ± 0.040 

7.51 ±0.035 

5.95 ± 0.027 

0.87 

8.27 ± 0.038 

7.13 ±0.033 

5.63 ± 0.026 

0.66 

7.33 ±0.034 

6.34 ± 0.029 

5.02 ± 0.023 

0.48 

6.54 ± 0.030 

5.62 ± 0.026 

4.36 ± 0.020 

0.35 

5.94 ± 0.027 

5.08 ± 0.023 

3.79 ±0.017 

0.25 

5.44 ± 0.025 

4.63 ±0.021 

3.40 ±0.016 

P 

Spectrum: Pereira et al. (2009b). EW: this work 

1.00 

8.63 ± 0.060 

7.55 ±0.052 

5.99 ±0.041 

0.816 

8.08 ± 0.056 

7.03 ± 0.049 

5.54 ±0.038 

0.608 

7.22 ± 0.050 

6.25 ± 0.043 

4.85 ± 0.034 

0.424 

6.33 ± 0.044 

5.46 ±0.038 

4.16 ±0.029 

0.197 

4.94 ± 0.034 

4.16 ±0.029 

3.02 ±0.021 

P 

Spectrum 

and EW: Pereira et al. (2009b) 

1.00 

8.48 ± 0.09 

7.44 ± 0.09 

5.90 ± 0.09 

0.816 

8.00 ±0.14 

6.96 ±0.15 

5.47 ±0.11 

0.608 

7.21 ± 0.09 

6.21 ±0.09 

4.83 ± 0.09 

0.424 

6.38 ±0.10 

5.44 ±0.10 

4.17 ±0.10 

0.197 

5.09 ±0.10 

4.22 ± 0.33 

3.08 ± 0.28 



Fig. 4. Equivalent widths as measured in this work from WCLC 
and Pereira’s spectra (red and black lines, respectively), com- 
pared to the measurements taken from the work of lPereira et ali 
(2009b) (black dots). See Table[2]for more details. 


4.2.1. Equivalent width measurements 

For a given model atmosphere, a grid of theoretical equivalent 
widths, WcaicCu), is computed as a function A(O) and Sh by in¬ 
tegrating the non-LTE synthetic line profile over a wavelength 
interval of +2 A, separately for each triplet component. Here the 
inclination angle g covers the values for both observed spec¬ 
tra (WCLC and Pereira). In addition, the equivalent widths of 
the disk-integrated line profiles were obtained by appropriate g- 
integration, denoted as VV _ ca ic in the following. 

The EWs of the observed line profiles, W^ s (g), W 0 b s , are ob¬ 
tained as a by-product of the 3D non-LTE line profile fitting (see 
Sect. 14. Il l: W 0 b s = Wcaic(A(0)*,5^), where A(O)* and 5^ are 
the best-fit parameters for the respective line profile, obtained 
by fitting each //-angle individually. This method of EW mea¬ 
surement ensures an optimal consistency between observed and 
theoretical equivalent widths. The errors of the observed EWs, 
cr 0 , were estimated from the formal errors of A(O) and .S'h for 
the best-fitting theoretical line profile. The EWs of the Oi IR 
triplet measured in this way are listed in Tab.[2] for all spec¬ 
tra; for the Pereira spectra, we also give the original measure¬ 
ments by IPereira etak ( 2009b ) (bottom). The EWs from both 
the WCLC and the Pereira spectra, as collected in Tab. [2] are 
plotted in Fig. [4] 

The EWs derived with our method from the Pereira spectra 
are generally slightly larger than those measured by the authors 
of the original work. This is because we integrate the synthetic 
single line profile over a very broad wavelength window. Since 
the Oi lines have extended Lorentzian wings, our EWs are sys¬ 
tematically higher. For the weakest lines at g — 0.197, the im¬ 
pact of the Lorentzian wings is diminished, and the two mea¬ 
surements are in very close agreement (see Table 0. We argue 
that, for the purpose of EW fitting, our EW measurement method 
is superior to direct integration of the observed line profiles be¬ 
cause it ensures consistency of theoretical and observed EWs. 

4.2.2. EW fitting of disk-center and full-disk spectra 

We first apply the method to determine S'h and A(O) from the si¬ 
multaneous fit of the disk-center and disk-integrated EWs in the 
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SH SH SH 

Fig. 5. Lines of constant equivalent width in the An - A(O) plane, as computed from the grid of 3D non-LTE synthetic line profiles 
for the three triplet components. The thick (red) solid line refers to the EW measured in the Neckel disk-center spectrum, while the 
thick (green) dashed line is for the measured EW in the Neckel disk-integrated (flux) spectrum. The solution (.S’n, A(O)) given in the 
legend is defined by the intersection point of the two lines (red diamond). The black diamonds define the uncertainty range of the 
solution, assuming an EW measurement error of w 0.5%, as indicated by the thin lines offset from the original .S'n - A(O) relations 
by ±AA(0) = 0.0025. 


Neckel intensity and flux spectra, respectively. For this purpose, 
we compute for each of the three triplet components two grids 
of equivalent widths, W(Ah, A(O); //= 1) and IV(.S'h, A(O); flux), 
from the corresponding sets of 3D non-LTE line profiles. These 
grids are then used to construct contour lines in the .S'h - A(O) 
plane where VV' ca ic = W 0 b s , as shown in Fig. [5] The thick (red) 
solid line and the thick (green) dashed line refers to the EW mea¬ 
sured in the Neckel disk-center and disk-integrated (flux) spec¬ 
trum, respectively (see Table 13. The slope of the contour line 
of the flux EW is somewhat steeper than that for disk-center, 
indicating that the lines in the flux spectrum are slightly more 
sensitive to changes in Ah- The unique solution (.S'h, A(O)) de¬ 
fined by the intersection point of the two lines (red diamond) is 
given in the legend of each panel. 

Since the angle at which the two lines intersect is rather 
small, the position of the intersection point is very sensitive to 
the exact values of the two equivalent widths. To quantify the 
uncertainties, we have also plotted in Fig. 0 (thin lines) the con¬ 
tour lines shifted up and down by AA(0)= +0.0025 dex, corre¬ 
sponding to errors in EW of less than 0.5%. The black diamonds 
show the solutions obtained by changing the two EWs by this 
small amount in opposite directions. With this definition of the 
error bars, we find A(O) = 8.77 ± 0.02, A H = 1.85 ± 0.65. 

These results are somewhat lower than the corresponding 
results obtained from line profile fitting (Sect. 14. 13l >. with a 
marginal overlap of the error bars. We conclude that the idea of 
simultaneously fitting intensity and flux spectra works in prin¬ 
ciple, but is highly sensitive to small errors in the measured 
EWs and to systematic errors in modeling the // dependence of 
the oxygen triplet lines. The results obtained with this method 
are therefore only considered a rough indication of the range 
of A(O) compatible with our modeling. As we see below, a 
more accurate determination of the oxygen abundance is pos¬ 
sible when fitting simultaneously the line profiles measured at 
different //-angles. 


4.2.3. EW fitting of the //-dependent WCLC and Pereira 
spectra 

We now determine Ah and A(O) from the simultaneous fit of the 
EWs at all //-angles, using the two independent sets of observed 
//-dependent spectra of WCLC and Pereira, respectively. As be¬ 
fore, we treat the three triplet components separately. In contrast 
to the method employed bv lPereira et al . ( 2009b ). we do not as¬ 
sign a special role to the fit at disk center. The idea of our method 
is to find the minimum^ 2 in the Ah - A(O) plane, where x 1 is 
defined as 


,T(Ah,A(0)) 


j=i 


(w° bs - W; alc (A H , A(O))) 2 


° 2 u +a ij 


(5) 


We sum over the different //-angles (J = 1 ... m); W /obs and W calc 
are the measured and computed equivalent width, respectively; 
cr y is the measurement error of U /ohs as given in Table|3 and <X 2 
is the observational error due to insufficient time averaging given 
the limited field of view covered by the slit on the solar surface. 
The error <x 2 is only considered in case of the WCLC spectra, 
where the exposure times are less than the typical granule life¬ 
time and the slit covers roughly three times the surface area of 
the 3D model. We derive <x 2 from a bootstrap-like test, where 
a random subsample of three different snapshots is repeatedly 
drawn from a larger sample of 20 3D snapshots without replace¬ 
ment. For each //-value, we generated 10 000 realizations of the 
possible appearance of the solar surface in the field of view, us¬ 
ing the same A(O) - Ah combination as for the equivalent width 
measurement (see above). The trial-to-trial variance of the de¬ 
rived distribution of the subsample averaged EWs is then taken 
as cr 2 . Pereira’s spectra have much better statistics (large num¬ 
ber of observations per //-angle), and there is no need to consider 
this kind of measurement error in the case of these spectra. 

For the purpose of EW fitting, we compute a grid of equiva¬ 
lent widths, W(Ah, A(O); //) from the corresponding sets of 3D 












12 


M. Steffen et al.: Solar oxygen abundance determined from the 111 triplet lines 



1.000 

4.000 

16.00 

64.00 

256.0 

1024 

4096 



1.000 

4.000 

16.00 

64.00 

256.0 

1024 

4096 


Fig. 6. Maps of^ 2 (color-coded) in the .S'h - A(O) plane, indicating the quality of the simultaneous fit over all //-angles of observed 
and computed 3D non-LTE EW, considering each of the three triplet components separately. The location of the minium x 2 is 
marked by a white plus sign, while the white ellipses bound the regions where x 2 ~ X^in < ' (solid) and 2.3 (dashed). The latter 
corresponds to the simultaneous confidence region of .S'n and A(O) containing 68.3% in the case of normally distributed data. 
Contour lines indicating where W’ ohs = W ca i c are superimposed for each //-angle. The top and bottom panels show the results for the 
WCLC and Pereira spectra, respectively. Line identification and best-fit parameters are given at the top left of each panel. 


non-LTE line profiles, where 0 < 5 h < 8/3, 8.65 < A(O) < 
8.83, and // covering the observed limb angles. This allows us 
to compute x 2 (Sn, A(O)) in the .S'h - A(O) plane, separately for 
each of the three triplet components, as shown in Fig.[6]for both 
observational data sets. The best-fit solution is defined by the lo¬ 
cation of the minimum ofx 2 , as indicated in the plots by a white 
plus sign, which is found by EW interpolation using third-order 
polynomials. 

The results of the EW fitting for various model atmospheres 
are summarized in Tabled listing A(O) and .S’h of the best fit 
for the WCLC and Pereira observed spectra. The formal fitting 
errors of A(O) and .S’n are given by the projection of the contour 
line defined by x 2 ~X 2 „ m = 1 (white ellipse in Fig. [6} on the A(O) 
and Sh axis, respectively. We also provide the reduced^ 2 of the 
best fit, which is defined as^ 2 ed = x 2 /( ,n ~ p)* where p - 2 is the 
number of free fitting parameters, as a measure of the quality of 
the fit. 

The high quality of the EW fit achieved with the 3D model is 
also illustrated in Fig. [7] where we show the center-to-limb varia¬ 
tion of observed and theoretical EWs for the best-fitting parame¬ 
ters and for the two extreme values Sh= 0 and 8/3. Corresponding 
plots for all ID model atmospheres investigated in this work are 
provided in Appendix[F| 

To derive the best estimate the solar oxygen abundance from 
EW fitting, we rely on the results of the 3D model atmosphere, 
which show the smallest line-to-line scatter and the lowest^ 2 . 


Computing the mean and errors for the two data sets as before in 
Sect.ECO we obtain <A(0)> = 8.754 ± 0.010dex (WCLC data 
set) and (A(O)) = 8.771 ± 0.006 dex (Pereira data set). Again, 
the formal 1 cr error bars of these two abundance estimates do 
not overlap (by a narrow margin). Nevertheless, we derive the 
best estimate by averaging over the two data sets, arriving at 
«A(0)» = 8.763+0.012, and «S H » = 1.6+0.2. These numbers 
are fully consistent with the final result obtained from the line 
profile fitting described in Section l4.1.21 As in the latter case, the 
best-fit solutions, derived from the (3D) and ID LHD modeling, 
turn out to be unphysical (negative .S'h). We discard these results 
as irrelevant and exclude them from Tabled (see also Sect. 15.21 
below). 


5. Discussion 

5.1. Comparison with the work of Pereira et al. (2009) 

In the present work, we simultaneously derived the oxygen abun¬ 
dance, A(O), and the scaling factor for collisions with neu¬ 
tral hydrogen, .S'h , using two different methods of analyzing 
the //-dependent spectra of the O i IR triplet: line profile fitting 
and equivalent width fitting. While the former method yields 
A(O) = 8.757 ± 0.010, S H = 1.6 ± 0.2, the latter gives A(O) = 
8.763 + 0.012, Sh = 1.6 + 0.2. In both methods, the spectra at 
the different //-angles are assigned equal weight. 
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Table 3. Results of non-LTE EW fitting with various model atmospheres: A(O), .S'h- values, and reduced^ 2 of the best fit for two 
different observed spectra (WCLC and Pereira 2009); the difference between A(O) derived from the two different spectra is given in 
Col. (4). Results from the (3D) and LHD model atmospheres are excluded because no minimum in^ 2 is found in the range .S'n > 0. 



A(O) 

AA(O) 


?H 

xL 


Data set: 

WCLC 

Pereira 


WCLC 

Pereira 

WCLC 

Pereira 

A, nm 

3D model atmosphere 

777.2 

8.763 ± 0.018 

8.765 ±0.014 

-0.002 

1.63 ±0.35 

1.67 ±0.25 

0.831 

0.500 

777.4 

8.743 ± 0.021 

8.772 ±0.013 

-0.029 

1.28 ±0.33 

1.69 ±0.25 

0.448 

0.457 

777.5 

8.757 ± 0.017 

8.776 ±0.010 

-0.019 

1.56 ±0.36 

1.96 ±0.26 

0.235 

0.401 

Mean 

8.754 ± 0.010 

8.771 ±0.006 

-0.017 

1.49 ±0.19 

1.77 ±0.17 

- 

- 

A , nm 


HM model atmosphere, = 

0.8 km/s 



777.2 

8.699 ± 0.030 

8.677 ±0.031 

0.022 

0.55 ± 0.25 

0.38 ±0.16 

2.135 

0.989 

777.4 

8.689 ±0.031 

8.705 ± 0.022 

-0.016 

0.42 ± 0.22 

0.49 ±0.17 

1.414 

0.477 

777.5 

8.731 ±0.022 

8.748 ±0.014 

-0.017 

0.73 ±0.26 

0.90 ±0.19 

0.371 

0.159 

Mean 

8.706 ± 0.022 

8.710 ±0.035 

-0.004 

0.57 ±0.16 

0.59 ± 0.28 

- 

- 

A, nm 


HM model atmosphere, £ m i C = 

1.2 km/s 



777.2 

8.670 ± 0.044 

8.650 ± 0.026 

0.020 

0.50 ±0.21 

0.36 ± 0.08 

2.330 

1.371 

777.4 

8.660 ± 0.030 

8.674 ± 0.026 

-0.014 

0.37 ±0.16 

0.42 ±0.16 

1.578 

0.695 

777.5 

8.705 ± 0.023 

8.721 ±0.015 

-0.016 

0.65 ± 0.25 

0.80 ±0.18 

0.420 

0.208 

Mean 

8.678 ± 0.024 

8.682 ± 0.036 

-0.004 

0.51 ±0.14 

0.52 ± 0.24 

- 

- 




i.i_i_i_i_i_i_i 

0.2 0.4 0.6 0.8 1.0 

P 


Fig. 7. Top panels: center-to-limb variation of observed (symbols) and theoretical EWs derived from the 3D model atmosphere 
for 5 h= 0 (blue dashed), 8/3 (red dotted) and best-fit value (black solid), shown individually for the three triplet components, with 
wavelength increasing from top to bottom. For each .S'n, the corresponding A(O) is defined as the value that minimizes x 2 (Eq.[5} 
at fixed Sh- Bottom panels: abundance difference A (p) - A, where At//) and A denote A(O) obtained from line profile fitting at 
individual //-angles and at all //-angles simultaneously. Left and right panels refer to the WCLC and Pereira spectra, respectively. 


IPereira et alJ ( 2009 b) followed a different methodology, us¬ 
ing a mix of EW fitting and line profile fitting, and assigning 
a distinguished role to the fits at disk center. More precisely, 
they first derive .S'n from simultaneously fitting the equivalent 
widths at all //-angles and all three triplet components. For a 
given .Sh, the oxygen abundance is varied to find the local min¬ 
imum / y 2 (.S'n). Then the minimum of this function defines the 
global minimum of % 2 , which corresponds to their best-fitting 


global value of .S'n = 0.85. Fixing .Sh at this value, these authors 
determine the oxygen abundance by fitting the line profiles at 
disk center, separately for each triplet component, and average 
the resulting A(O) values to obtain their best estimate of the oxy¬ 
gen abundance, A(O) = 8.68. 

In Fig. [8] we show a comparison of the A(O) - .S'n rela¬ 
tion obtained from our 3D-NLTE modeling with the results of 
[Pereira et al.l ( 2009bl) . for each of the three triplet lines and the 
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two different disk-center spectra. For a given Sh, the oxygen 
abundance is determined from disk-center line profile fitting, as 
described in Sect. l4.1.T1 (first sweep, but with fixed Sh). The fig¬ 
ure immediately shows that_our A(O) - .S'n relation differs sig¬ 
nificantly from that of iPereira et al.l (l2009bl . their Table 3), im¬ 
plying a systematically higher A(O) for given .S'n. For the value 
,S[| = 0.85, which was determined bv IPereira et al.l (1200% ). our 
relation suggests A(O) =» 8.71 - 8.72. This value is clearly in¬ 
compatible with the value A(O) « 8.64-8.68, found by the latter 
authors. However, we are cautious about making a direct com¬ 
parison of the A(O) - Sh relations since possible differences in 
the particular implementation of the Drawin formula (definition 
of Sh) might introduce a bias (see Sect. l5.3.TI >. This is the case 
with the exception of Sh = 0 and oo (LTE). 

In addition, the line profile fits that we achieve with our 3D 
non-LTE synthetic line profiles (see Fig.|3} are superior to the fit¬ 
ting results shown in the work of lPereira et alj d2009bl) (see their 
Fig. 6 and Sect. 4.1.2.) where, generally, a less perfect agreement 
was found. 

We can only speculate about what causes the syst ematic dif¬ 
ference between our results and those obtained bv iPereira et al.l 
( 2009bi) . Differences in the employed oxygen model atom, the 
adopted atomic data (like photoionization and collisional cross 
sections), and other numerical details entering the solution of the 
non-LTE rate equations, are among the prime suspects. Since we 
even see systematic A(O) differences of about 0.03 dex in LTE 
(right column of Fig. [8), we conclude that differences in the 3D 
hydrodynamical model atmospheres and / or the spectrum syn¬ 
thesis codes must also play a role. 

It seems very likely that the above mentioned differences in 
the theoretical 3D plus non-LTE modeling are the main rea¬ 
son for the discrepancy in the final oxygen abundance of al¬ 
most/). 1 dex between the present work and that of lPereira et al.l 
(l2009bl) . The different approaches of analyzing the data consti¬ 
tute an additional source of uncertainty. 

5.1.1. Impact of equivalent width measurements 

As described in Sect. l4.2.Tl we systematically derive larger 
equivalent widths than IPereira et al.l ( 2009b ) near the disk center 
and lower EWs near the limb, from the same observed spectra 
(see Table[2} due to different measurement methods. One might 
therefore ask how much of the difference between the oxygen 
abundance derived in the present work and by Pereira et al. can 
be ascribed to differences in the measured equivalent widths. To 
answer this question, we have repeated our 3D, non-LTE, EW 
fitting procedure with the original equivalent widths as given 
bv IPereira et al . (2009b ). We find that the derived oxygen abun¬ 
dance decreases by AA(O) » -0.045 dex, explaining roughly 
half the discussed A(O) discrepancy. We point out, however, that 
the almost perfect consistency between profile fitting and EW fit¬ 
ting results is destroyed if we were to adopt the equivalent widths 
measured bv IPereira et al.l (2009b) . Note also that the mismatch 
illustrated in Fig.[8]is completely independent of any EW mea¬ 
surements. 

5.2. Interpretation of the ID results 

As mentioned in Sect. 14. 11 the oxygen abundance derived from 
line profile fitting based on ID models (in particular (3D) and 
LHD) is significantly lower, by up to -0.15 dex, than the value 
derived with the 3D modeling. A very similar discrepancy be¬ 
tween ID and 3D results is obtained from the method of equiva- 



Fig. 8. Oxygen abundance determined from disk-center line pro¬ 
file fitting as a function of .S'h for the three triplet lines and 
the two different observed spectra. The 3D-NLTE results of the 
present work (solid lines, symbols) are compared with those of 
IPereira et al.l ( 2009b . Tab. 3, dashed lines), obtained from their 
profile fits of the same disk-center spectra. 



Fig. 9. Minimum x 2 as a function of Sh obtained from global 
EW fitting with the (3D) and LHD model atmospheres. For each 
of the three triplet lines, we plot the minimum of x 2 (nonre- 
duced) at fixed Sh- No minimum of x 2 is found in the range 
0 < Sh < 8/3. At Sh=0, the implied oxygen abundances are 
A(O) = 8.61, 8.62, 8.63 «3D)), and A(O) = 8.59, 8.60, 8.61 
(LHD), for A 111 .2, 111 A, and 777.5 nm, respectively. The best 
fits would be obtained at negative Sh and even lower A(O). 


lent width fitting. The main reason is that the ID models predict 
a steeper gradient of the line strength with respect to limb angle 
p (compare Fig. [7] with Figs. IF.il and IF.3I >. Under such circum¬ 
stances, the best global fit of the observed center-to-limb varia¬ 
tion of the triplet lines would formally require a negative Sh, as 
demonstrated in Fig. [9] Because of the tight correlation between 
Sh and A(O) (see Fig. [8}, a low Sh implies a low oxygen abun¬ 
dance: while Sh ~ 1.6 for the 3D model, Sh = 0 for the (3D) and 
LHD models, explaining the exceedingly low oxygen abundance 
derived with these lD_models. We add that our findings are in 
line with the results of lPereira et al.l ( 2009bl their Fig. 5). 
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5.2.1. The microturbulence problem 

Since the triplet lines are partially saturated, their formation is 
sensitive to the choice of the microturbulence parameter, £ m ; c , 
to be used with the ID models. While the value of £ m j c can 
be determined empirically, the main issue in the present con¬ 
text is the common assumption that 4 ; rn i C is independent of //. In 
fact, it is well known empirically that both micro- and macro¬ 
turbulence increase towards the solar limb (e.g. jHolweger et all 
197K). Evidently, we could improve the quality of the ID fits 
by introducing a //-dependent microturbulence parameter. With 
a more realistic £- m ic increasing towards the limb, the slope of 
W(jj) at fixed A(O) would diminish, such that the best fit would 
require a larger (positive) .S'h and hence a higher oxygen abun¬ 
dance. Although it may be worthwhile, we have refrained from 
exploring the effects of a //-dependent microturbulence on the 
results of the ID models in the present work. 

It is obvious that the ID models with a fixed microturbu¬ 
lence are unrealistic, and lead to unphysical conclusions about 
the scaling factor (S'h < 0). The oxygen abundance derived from 
these ID models for the Oi IR triplet lines is therefore system¬ 
atically biased towards exceedingly low values. Even though the 
analysis based on the HM model atmosphere yields positive .S’h 
values, the resulting A(O) is considered as misleading, too, since 
the HM model suffers from the same microturbulence issues as 
the (3D) and LHD models, albeit in a less obvious way, because 
the problems are somehow masked by the empirical temperature 
structure of the HM model. Consequently, we have ignored all 
ID results for the estimation of the final solar oxygen abundance. 
We rather rely on the analysis with the 3D model, which has a 
built-in hydrodynamic al velocity field that properly represents a 
//-dependent micro- and macroturbulence. 

5.2.2. 1D - 3D comparison 

As described above, the microturbulence problem renders a di¬ 
rect comparison of the abundance determinations from ID and 
3D models meaningless. Nevertheless, it is possible to ask what 
oxygen abundance would be obtained from the ID models if we 
use the same .S'h- value as found in the 3D analysis, i.e., consider¬ 
ing this .S'n-value as a property of the model atom. Fixing .S'h in 
this way, we determined A(O) from the ID models by fitting the 
observed equivalent width at disk center, assuming the canonical 
microturbulence value of £ m j c = l.Okm/s. The results are sum¬ 
marized in Table[4] 

Remarkably, the HM model atmosphere requires almost ex¬ 
actly the same oxygen abundance as the 3D model to match the 
measured disk-center equivalent width. This is indeed what is 
expected from a semi-empirical model. The averaged 3D model 
gives only a slightly lower A(O), suggesting that the effect of 
horizontal inhomogeneities is of minor importance, at least for 
the line formation at // = 1. In contrast, the LHD model indicates 
a significantly lower oxygen abundance, which we attribute to 
the somewhat steeper temperature gradient in the line forming 
region (see Sect. l3.1.3l >. 

5.3. Treatment of inelastic coiiisions by neutral hydrogen 

Collisions of neutral oxygen atoms with neutral hydrogen give 
rise to bound-bound (excitation) and bound-free (ionization) 
transitions. As shown previously, they play a very important role 
in establishing the statistical equilibrium level population of O i. 
Unfortunately, quantum mechanical calculations or experimen¬ 
tal measurements of the relevant collisional cross sections are 


Table 4. Results of non-LTE EW fitting with various model at¬ 
mospheres for given disk-center equivalent width W\ (as mea¬ 
sured from the Pereira data set, see Table[2} and .S'h- factor (as 
obtained from the global EW fit with the 3D model, see Table[3}. 
The oxygen abundance A(O) required to reproduce W\ with the 
3D, averaged 3D, HM, and LHD model, respectively, is listed in 
Cols. (4) - (7). For all ID models, a microturbulence of /- m j c = 
l.Okm/s is assumed. 


A 

[nm] 

Wi 

[pm] 

Sh 

3D 

A(O) 

(3D) HM 

LHD 

777.2 

8.63 

1.67 

8.758 

8.744 

8.753 

8.686 

777.4 

7.55 

1.69 

8.766 

8.753 

8.765 

8.695 

777.5 

5.99 

1.96 

8.770 

8.759 

8.778 

8.703 

Mean 

- 

- 

8.765 

8.752 

8.765 

8.695 


so far unavailable. We have therefore resorted to calculating the 
excitation and ionization rates due to collisions by neutral hydro¬ 
gen by the classical Drawin formula. Following common prac¬ 
tice, we scaled the resulting rates by an adjustable factor, .Sh, 
which is determined from the best fit to the observations together 
with the oxygen abundance, as detailed above. 

5.3.1. Alternative implementations of the Drawin formula 

In the present work, the cross sections for inelastic collisions 
with neutral hydrogen are computed according to the classical 
recipe of Drawin (lDrawinlll969l: fSteenboc k & Holwei/erl ll984; 
Lambert 1993), and hence are proportional to the oscillator 
strength / for all bound-bound transitions. For optically forbid¬ 
den transitions, the collisional cross section is therefore zero. 
All hydrogen related collision rates included in our oxygen 
model atom, bound-bound and bound-free, are then scaled by 
the sam e factor .S’n. It is not clear whether other authors (e.g., 
IPereira et al.ll2009b ) might have adopted different assumptions. 
We have therefore investigated two test cases to understand the 
consequences for the derived oxygen abundance of (i) allowing 
for collisional coupling through optically forbidden transitions, 
and (ii) assuming a constant collision strength independent of / 
for all allowed transitions. 

To reduce the computational cost, we have restricted the test 
calculations to EW fitting of the Pereira data set, using the HM 
ID model atmosphere with £ m j c = 1.2 km/s. The results are sum¬ 
marized in Table[5] The oxygen abundance and .S’h- factor ob¬ 
tained with the standard setup are shown in Cols. (2) and (3) 
(‘Case O’, cf. Tabled. 

In ‘Case 1’, we have included collisional excitation for all 
possible transitions according to the Drawin formula, where the 
cross sections of the optically allowed transitions are propor¬ 
tional to the oscillator strength / (as in ‘Case O’). In this case, 
an effective minimum oscillator strength of / = 0.001 is adopted 
for all forbidden transitions, including the transitions between 
the triplet and quintuplet system, as suggested by iFabbian et al.l 
d20m Scaling by S H is restricted to the allowed transitions. 
The effect of the additional collisional coupling via the forbidden 
transitions on the oxygen abundance derived from EW fitting of 
the triplet lines is marginal: on average, A(O) is reduced by only 
0.004 dex. 

In ‘Case 2’, the treatment of the collisional coupling in the 
forbidden transitions is the same as in ‘Case 1 ’, but the cross sec¬ 
tions of all optically allowed transitions are computed with f — 1. 
Again, scaling by S’h is restricted to the allowed transitions. The 
effect of removing the proportionality between collision strength 
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Table 5. Resulting oxygen abundance A(O) and scaling factor 
.S '11 derived from non-LTE EW fitting of the Pereira data set, us¬ 
ing the HM model atmosphere with i : rn i c = 1.2 km/s. Cases (0) - 
(2) refer to different assumptions about the excitation rates due 
to collisions by neutral hydrogen (see text for details). 


A 

[nm] 

Case 0 
A(O) S‘ 0) 

Case 1 

AA(O) 

Case 2 

AA(O) 

777.2 

8.650 

0.36 

-0.004 

0.818 

-0.011 

0.515 

777.4 

8.674 

0.42 

-0.005 

0.854 

-0.014 

0.512 

777.5 

8.721 

0.80 

-0.003 

1.000 

-0.006 

0.722 

Mean 

8.682 

0.52 

-0.004 

0.891 

-0.010 

0.583 


and oscillator strength for the allowed transitions on the derived 
oxygen abundance amounts to « -0.01 dex, while the numerical 
value of the required .S'h is reduced by roughly a factor 2 with 
respect to ‘Case O’. 

The effects shown in Table[5] are expected to be of similar 
magnitude in the full 3D analysis. We conclude from this ex¬ 
periment that the oxygen abundance derived from the non-LTE 
analysis of the O i IR triplet lines is not particularly sensitive to 
the details of the implementation of the Drawin formula, even if 
the numerical value of .S'h can differ significantly. 


5.3.2. Assumption of a universal S H 

Although this approach is commonly accepted, the assumption 
of a universal scaling factor is questionable. It could well be that 
each line is characterized by a slightly different value of .S'h, and 
/ or that individual bound-bound and bound-free transitions re¬ 
quire significantly different scaling factors. This would not be 
surprising, as it is well known that the classical treatment of col¬ 
lisions does not have a rigorous physical foundation. In fact, such 
nonuniform scaling has recently been demonstrated for Mg + 
H by Barkle m et al.l (1 201 2l), In princ iple. .S'u should depend on 
temperature (see Barklem et al. 2011) and, because the different 
triplet lines form at different atmospheric heights, S’h could be 
expected to vary slightly from line to line. Our analysis allows 
for this possibility, since we treat each of the triplet components 
separately. However, we do not find any systematic trend in .S'h 
with wavelength, and consider the line-to-line variation of .S'h 
as random, reflecting the intrinsic uncertainties of our method. 
With this reasoning, a systematic variation of .S'h as a function 
of the limb angle also seems unlikely, and our analysis assumes 
a //-independent .S'h- In summary, the treatment of the collisions 
by neutral hydrogen in terms of the empirical scaling factor .S'h 
appears unsatisfactory. 


6. Conclusions 

The present analysis of the O i IR triplet is based on two data sets 
of intensity spectra recorded at different /./-angles across the solar 
disk: a set of spectrograms taken by one of the authors (WCL) 
with the rapid scan double-pass spectrometer at the McMath- 
Pierce Solar Telescope, Kitt Peak, on_12 September 2006, and 
the set of high-quality spectra used by Pereira et al.l (l2009bt) . ob¬ 
served with the Swedish 1-m Solar Telescope on Roque de Los 
Muchachos, La Palma, in May 2007. 

The photospheric oxygen abundance, A(O), and the scaling 
factor for collisional excitation / ionization by neutral hydro¬ 
gen, .S’h, were determined simultaneously, using two different^- 2 
fitting methods of minimizing the difference between observed 
and synthetic spectra. In the first method, we simultaneously 


fit the line profiles at all /./-angles. Averaging the results of the 
three triplet components and the two different sets of observed 
spectra, we find an oxygen abundance of A(O) = 8.757 + 0.010 
(Sh = 1.6 ± 0.2), where the error is defined by the standard de¬ 
viation of the individual results. In the second method, we fit 
equivalent widths instead of line profiles, and obtain a perfectly 
consistent result, A(O) = 8.763 ± 0.012 (Sh = 1.6 ± 0.2). 

The above numbers were derived by employing 3D non-LTE 
synthetic spectra, consistently taking non-LTE line formation in 
a 3D hydrodynamical model of the solar atmosphere into ac¬ 
count. When performing the same analysis with synthetic spec¬ 
tra computed from various ID model atmospheres, we obtain 
much less perfect line profile fits, a larger line-to-line variation 
of A(O) and S'h, and a significantly lower value of the best-fitting 
oxygen abundance. We argue that the ID results are systemati¬ 
cally biased towards low Sh and A(O) values, mainly as a con¬ 
sequence of the assumption of a /./-independent microturbulence. 
Consequently, we have ignored all ID results for the estima¬ 
tion of the final solar oxygen abundance in the present work. 
We rather rely on the analysis with the 3D model, which has a 
built-in hydrodynamical velocity field that properly represents a 
//-dependent micro- and macroturbulence. 

All individual results derived from the various observations 
of the O i IR triplet lines with different analysis methods, fall in 
the range 8.74 < A(O) < 8.78. We thus arrive at a final 3D non- 
LTE oxygen abundance of A(O) = 8.76 ± 0.02. This coincides 
with the value recommended in Paper I (see also ICaffau et al.l 
1201 li) based on various atomic Oi lines, but is significantly 
higher than suggested by the analysis of iPereira et alJ (l2009bh . 
who find A(O) = 8.68. The reason for this discrepancy of al¬ 
most 0.1 dex is unknown. Presumably, differences in the 3D hy¬ 
drodynamical model atmospheres (C05B0LD versus STAGGER), 
in the employed oxygen model atom, or in the adopted atomic 
data may offer an explanation. 

The present detailed investigation leads us to conclude that, 
from the CO 5 BOLD perspective, a solar photospheric oxygen 
abundance lower than A(O) = 8.70 now appears very unlikely. 
In view of the recent laboratory measurements of iBailev et al.l 
(12015h . indicating that the iron opacity near the solar radiation / 
convection zone boundary is significantly higher than assumed 
so far, the conflict between surface chemical abundances and 
helioseismic models might thus be less severe than previously 
thought. 

Finally, we note that the treatment of the collisions by neu¬ 
tral hydrogen in terms of an empirical scaling factor .S'h is far 
from satisfactory. The assumption that the same .S’u is valid for 
all transitions is questionable and introduces an additional sys¬ 
tematic uncertainty of unknown magnitude in the determination 
of the photospheric solar oxygen abundance from the Oi IR 
triplet. Clearly, accurate theoretical or experimental collisional 
cross sections for neutral oxygen are badly needed. 
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Appendix A: The oxygen model atom designed for 
solar applications 

In Tables IA. II and IA.2I we provide details about the 
22 energy levels and 54 radiative transitions consid¬ 


ered in the oxygen model atom used in this work. The 
data are taken from the NIST Atomic Spectra Database 
(http: //www. nist. gov/pral/data/asd. cfm/l. 

Table A.l. Atomic energy levels of the oxygen model atom used 
in this work. 


(1) 

Level 

# 

(2) 

Configuration 
and term 

(3) (4) 

Energy 

[Ryd] “ [eV] 

(5) 

Statistical 

weight 

1 

2p 3 P 

0.0000000 

0.0000 

9 

2 

2p ‘D 

0.1438880 

1.9577 

5 

3 

2p‘S 

0.3079407 

4.1898 

1 

4 

3s 5 S° 

0.6722253 

9.1461 

5 

5 

3s 3 S" 

0.6990967 

9.5117 

3 

6 

3p 5 P 

0.7894125 

10.741 

15 

7 

3p 3 P 

0.8069543 

10.979 

9 

8 

4s 5 S” 

0.8700480 

11.838 

5 

9 

4s 3 S° 

0.8768672 

11.930 

3 

10 

3d 5 D° 

0.8877632 

12.079 

25 

11 

3d 3 D° 

0.8883803 

12.087 

15 

12 

4p 5 P 

0.9030074 

12.286 

15 

13 

4p 3 P 

0.9076500 

12.349 

9 

14 

3s 3 D'’ 

0.9216869 

12.540 

15 

15 

5s 5 S° 

0.9305559 

12.661 

5 

16 

5s 3 S” 

0.9332467 

12.698 

3 

17 

4d 5 D° 

0.9373803 

12.754 

25 

18 

4d 3 D° 

0.9377701 

12.759 

15 

19 

4f 5 F 

0.9383157 

12.766 

35 

20 

4f 3 F 

0.9383165 

12.766 

21 

21 

5p 5 P 

0.9443120 

12.848 

15 

22 

5p 3 P 

0.9465332 

12.878 

9 


Appendix B: Ignoring the fine-structure splitting of 
the triplet 

In Sect. l3.2.n we mentioned that the upper level of the Oi IR 
triplet (3p 5 P) is treated as a single superlevel in our oxygen 
model atom, thus ignoring the fine structure splitting of this term. 
Since the energy differences between fine-structure levels are 
negligible, the LTE population of this superlevel equals the sum 
of those of the individual sublevels. If the fine-structure levels 
are merged into one superlevel, the three transitions making up 
the triplet are merged into a single combined transition, and the 
£//-value of this combined transition is the sum of the ^/-values 
of the individual fine-structure components. 

The radiative transition probability (Einstein coefficient or, 
equivalently, gf- value) of a bound-bound transition plays a two¬ 
fold role in the non-LTE problem: firstly, the photoexcitation rate 
at given radiation field scales with gf\ and secondly, the line 
opacity is proportional to gf, and hence the radiation field of the 
transition is a (non-linear) function of the g/-value. Since the 
three triplet components are well separated in wavelength and 
partially saturated in the Sun, merging them leads to a strongly 
saturated single line where the mean radiation field in the line, 
/, is diminished below the value in any of the individual original 
triplet components. The photoexcitation rates in the combined 
transition would therefore be significantly underestimated. 

To avoid this systematic bias, we reduce the line opacity of 
the combined triplet transition by a factor 3. This is equivalent 
to using, only for the calculation of /, the g/-value of the inter¬ 
mediate triplet component instead of the sum of the ^/-values of 
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Table A.2. Radiative transitions included in the oxygen model 
atom used in this work. 


1) 

Line 

# 

(2) (3) 
Levels 
lo up 

(4) 

A 

[A] 

(5) 

Oscillator 
strength, / 

(6) (7) 

Transition probability 

A ji Bjj 

1 

1 

2 

6372.0 

3.57E-11 

1.07E-2 

3.80E+00 

2 

1 

5 

1303.5 

5.19E-02 

6.12E+8 

1.13E+09 

3 

1 

9 

1040.1 

9.17E-03 

1.70E+8 

1.60E+08 

4 

1 

10 

1026.5 

8.33E-13 

1.90E-3 

1.43E-02 

5 

1 

11 

1026.6 

2.00E-02 

7.63E+7 

3.45E+08 

6 

1 

16 

977.2 

3.30E-03 

6.93E+7 

5.41E+07 

7 

1 

17 

974.1 

5.62E-05 

1.43E+5 

9.18E+05 

8 

1 

18 

972.5 

1.37E-02 

5.84E+7 

2.24E+08 

9 

2 

3 

5578.9 

1.16E-09 

1.26E+0 

1.08E+02 

10 

2 

4 

1727.1 

7.76E-10 

1.74E+0 

2.24E+01 

11 

2 

5 

1641.3 

4.43E-07 

1.83E+3 

1.22E+04 

12 

2 

14 

1172.6 

8.51E-06 

1.38E+4 

1.67E+05 

13 

3 

5 

2324.8 

1.12E-08 

4.61E+0 

4.40E+02 

14 

3 

14 

1484.5 

8.97E-06 

1.81E+3 

2.23E+05 

15 

4 

6 

7775.5 

1.00E+00 

3.69E+7 

1.31E+11 

16 

4 

12 

3948.5 

3.42E-03 

4.89E+5 

2.27E+08 

17 

4 

21 

3348.2 

3.02E-04 

5.99E+4 

1.69E+07 

18 

5 

7 

8448.8 

1.03E+00 

3.22E+7 

1.46E+11 

19 

5 

13 

4369.5 

6.50E-03 

7.58E+5 

4.77E+08 

20 

5 

22 

3692.4 

5.67E-04 

9.30E+4 

3.50E+07 

21 

6 

8 

11302.5 

1.70E-01 

2.67E+7 

3.23E+10 

22 

6 

10 

9266.4 

9.54E-01 

4.45E+7 

1.48E+11 

23 

6 

15 

6456.8 

1.71E-02 

8.25E+6 

1.86E+09 

24 

6 

17 

6159.0 

7.22E-02 

7.62E+6 

7.46E+09 

25 

7 

9 

13168.2 

1.81E-01 

2.14E+7 

3.97E+10 

26 

7 

11 

11289.9 

9.67E-01 

3.09E+7 

1.81E+11 

27 

7 

14 

7992.0 

8.84E-04 

5.61E+4 

1.17E+08 

28 

7 

16 

7256.4 

1.74E-02 

6.72E+6 

2.11E+09 

29 

7 

18 

7004.1 

4.28E-02 

3.53E+6 

5.00E+09 

30 

8 

12 

27644.6 

1.47E+00 

4.29E+6 

6.84E+11 

31 

8 

21 

12267.7 

1.67E-02 

2.47E+5 

3.44E+09 

32 

9 

13 

28935.2 

1.61E+00 

4.10E+6 

8.03E+11 

33 

9 

22 

13076.9 

2.17E-02 

2.83E+5 

4.78E+09 

34 

10 

12 

59762.0 

1.59E-01 

4.97E+5 

1.60E+11 

35 

10 

19 

18026.1 

1.00E+00 

1.48E+7 

3.05E+11 

36 

10 

20 

18021.2 

1.96E-02 

4.79E+5 

5.93E+09 

37 

10 

21 

16114.7 

9.43E-04 

4.04E+4 

2.55E+08 

38 

11 

13 

45608.5 

1.98E-01 

9.86E+5 

1.57E+11 

39 

11 

19 

18244.1 

2.98E-02 

2.56E+5 

9.13E+09 

40 

11 

20 

18248.6 

1.02E+00 

1.47E+7 

3.14E+11 

41 

11 

22 

15670.2 

4.04E-03 

1.83E+5 

1.06E+09 

42 

12 

15 

33083.7 

3.01E-01 

5.51E+6 

1.67E+11 

43 

12 

17 

26514.7 

1.13E+00 

6.45E+6 

5.04E+11 

44 

13 

14 

68287.3 

1.43E-02 

1.36E+4 

1.56E+10 

45 

13 

16 

36617.4 

3.05E-01 

4.82E+6 

1.82E+11 

46 

13 

18 

30985.0 

1.15E+00 

5.04E+6 

5.85E+11 

47 

14 

20 

54797.7 

9.83E-03 

1.56E+4 

9.04E+09 

48 

14 

22 

36687.2 

1.98E-03 

1.64E+4 

1.22E+09 

49 

15 

21 

66237.7 

1.88E+00 

9.57E+5 

2.10E+12 

50 

16 

22 

68568.3 

1.94E+00 

9.17E+5 

2.23E+12 

51 

17 

19 

974107.5 

4.78E-02 

2.40E+2 

7.81E+11 

52 

17 

21 

131456.0 

3.14E-01 

2.02E+5 

6.93E+11 
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all components. For the calculation of the photoexcitation rates 
at given radiation field, however, the total p/-value (i.e., the sum 
of the c//-values of the individual fine structure components) is 
applied. This procedure guarantees both J and photoexcitation 
rates to be consistent with the real situation. 

In the following, we show that the adopted approach is equiv¬ 
alent to having a detailed multiplet transition with three lines, 
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Fig. B.l. Departure coefficients of the Oi IR triplet levels as a 
function of the optical depth for three cases of the fine-structure 
treatment (see text for details). The black solid line shows the 
results obtained when 3p 5 P is treated as a superlevel (case 1); 
colored lines show the results when fine-structure splitting is in¬ 
cluded; dashed lines show case 2 (no collisional coupling be¬ 
tween fine-structure levels); solid lines show case 3 (strong col¬ 
lisional coupling between fine-structure levels). Top and bottom 
panels refer to terms 3p 5 P (upper level of the triplet) and 3s s S° 
(lower level of the triplet), respectively. 

with highly efficient collisions between the fine-structure lev¬ 
els. The latter approach was utilized in numerous works (e.g., 
lKiselmanlll993b IPrzvbilla et alJl200Ct lAsnlund et al.ll2004t) and 
is assumed to be a valid representation of the real situation. To 
demonstrate the equivalence of the two methods, we constructed 
a model atom that includes the fine-structure splitting of the 3p 5 P 
level. We used a solar LHD model atmosphere to compute the de¬ 
parture coefficients for three cases: 

1. Using the standard model atom where 3p 5 P is treated as a 
superlevel and gf is reduced by 1/3 for the calculation of the 
line opacity; 

2. Using a model atom that includes the splitting of the 3p 5 P 
term and no collisional coupling between the fine-structure 
levels; 

3. Using a model atom that includes the splitting of the 
3p 5 P term and strong collisional coupling between the fine- 
structure levels. 

In the last two cases, no scaling of the line opacity is ap¬ 
plied. Except for the treatment of the 3p 5 P term and the transi¬ 
tions that connect to it, all three model atoms are identical (see 
Sect. 13. 2. Ti l. For the test calculations, the oxygen abundance was 
set to A(O) = 8.77 and S H = 1.0. 
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The results of this experiment are shown in Fig. IB. II Firstly, 
it is clear that the treatment of the level 3p 3 P does not signif¬ 
icantly influence the lower level population of the triplet. This 
means that non-LTE line opacity of the O i IR triplet is not sensi¬ 
tive to the details regarding the fine-structure of 3p 5 P. Secondly, 
it is evident that, when no collisional coupling is applied between 
fine-structure levels, they exhibit slightly different departure co¬ 
efficients. Our analysis (not detailed here) has shown that this 
is caused by the different line strengths of the three triplet com¬ 
ponents, which leads to different photoexcitation rates in each 
of the triplet transitions and, in turn, to different departure co¬ 
efficients. Finally, the most important result is that, when ef¬ 
ficient collisions between the fine-structure levels are applied, 
the departure coefficients of the individual sublevels thermalize 
to the same values that were computed with the superlevel ap¬ 
proach. We consider this as a validation of our simplified ap¬ 
proach, which reduces the number of lines for which the radi¬ 
ation transfer has to be computed and thus facilitates a faster 
solution (by about 60 %) of the non-LTE statistical equilibrium 
problem. 

Appendix C: Effect of horizontal resolution on the 
3D-NLTE synthetic spectrum 

In Sections f3.2.2l and l3.2.3l we have mentioned that to speed up 
the calculations, we applied a horizontal reduction of the model 
atmosphere, using only every third vertical column in each hor¬ 
izontal direction (i.e., discarding 90% of the data defining the 
3D model atmosphere). In this section, we briefly validate this 
approach. 

To estimate the effect of the reduced horizontal resolution, 
we performed non-LTE calculations (Sect. l3.22B and subsequent 
spectrum synthesis (Sect. 13.231) with a fully resolved single 3D 
model snapshot. In these calculations, we set A(O) = 8.77 and 
5 h = 1 -0. To see how the effect changes across the solar disk, 
spectrum synthesis was performed at two inclination angles: p = 
1.0 and p = 0.2. 

Fig- EH shows the results of this experiment. It is seen that 
the effect of the reduced horizontal resolution does not exceed 
0.1 % of the continuum intensity anywhere in the line profile. The 
effect on equivalent widths is even smaller, 6W/W ~ 7 x 10 4 
and 12 X 10~ 4 for p = 0.2 and p = 1.0, respectively. Hence we 
may conclude that the reduced horizontal resolution that we have 
imposed on the spectrum synthesis does not influence the results 
in any significant way. 

Appendix D: Influence of the mixing-length 
parameter on 1D-NLTE spectrum synthesis of 
the Oi IR triplet lines 

The O i IR triplet lines share a lower level of high excitation po¬ 
tential and, hence, are expected to form in the deep photosphere. 
It could therefore be expected that 1D-NLTE spectral line syn¬ 
thesis of such lines is influenced by the choice of the mixing- 
length parameter, cumlt* which quantifies the convective effi¬ 
ciency and thus controls the temperature structure of the lower 
photosphere. To probe the possible influence of omlt, we com¬ 
puted a grid of 1D-NLTE line profiles using solar LHD model 
atmospheres computed with the different values of crMur- This 
grid of synthetic spectra is used to perform the same line profile 
fitting as described in Sect. 14.11 

Despite the fact that the fits are quite bad in general (see 
discussion in Sect. 14.11 and Fig. IE. 31) . we arrived at very simi- 
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Fig.C.l. Effect of reduced horizontal resolution on the 3D- 
NLTE spectral line profile of the middle component of the O i IR 
triplet. Top panels : synthetic spectral lines computed with full 
resolution (black line) and with resolution reduced by a factor 3 
in each horizontal direction (red dots). Bottom panels: difference 
in residual intensity between synthetic spectral lines computed 
with full and reduced horizontal resolution. A negative differ¬ 
ence means that the spectral line computed with full resolution 
is stronger at the particular wavelength. Left and right panels 
show the results at p - 0.2 and p - 1.0, respectively. 


lar results regardless of the chosen ctmlt- Table IdTTI shows the 
results of line profile fitting with different values of (Tmlt and 
fixed = 1.2km/s (assuming £ m i c = 0.8km/s yields simi¬ 
lar conclusions). Varying ckmlt between 0.5 and 1.5 changes the 
best-fitting A(O) by at most 0.006 dex (for the weakest triplet 
component). Also, the reduced^ 2 values of the best fit are only 
marginally influenced. We conclude that A(O) is very weakly 
dependent on the choice of ctmlt, and the goodness of the fits 
cannot be improved by modifying o:mlt- Hence our results de¬ 
rived from LHD models with otmlt = 10 should be represen¬ 
tative for all theoretical ID mixing-length models of the solar 
atmosphere. 


Table D.l. Best-fitting A(O) derived from the two sets of ob¬ 
served spectra, using ID LHD model atmospheres with different 
ctmlt values between 0.5 and 1.5, and fixed £' mlc = 1.2km/s. 
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Appendix E: Supplementary plots: line profile 
fitting results for ID model atmospheres 

In this appendix, we show the observed center-to limb variation 
of the line profiles of the O i IR triplet, separately for each com¬ 
ponent, in comparison with the best-fit theoretical non-LTE line 
profiles derived from various ID model atmospheres. The best- 
fit values of the oxygen abundance, A(O), are given in the legend 
of each plot, as well as in TableQ] (together with the best-fit .S'h 
values). The best fit correspond to the global minimum of as 
defined by Eq. (|4}, i.e., simultaneously fitting the line profiles at 
all //-angles. The quality of the fit can be judged from the re¬ 
duced x 2 values provided in the Table, as well as from the plots 
of the intensity difference 7 obs - 7 calc displayed below the plot of 
the line profiles. 


Appendix F: Supplementary plots: equivalent width 
fitting results 

In this appendix, we show the center-to-limb variation of the 
observed equivalent widths, separately for each component of 
the O i IR triplet, in comparison with the theoretical results de¬ 
rived from the 3D and various ID model atmospheres. For given 
.S’h= 0 and 8/3, the oxygen abundance is computed as the A(O) 
value that minimizes the global^/ 2 given by Eq. <[5}. The best fit 
is defined as the .S'h - A(O) combination that corresponds to the 
global minimum of x 1 - The best-fit values of Sh and A(O) are 
given in Tabled The quality of the fit can be judged from the re¬ 
duced x 2 values provided in the Table, as well as from the plots 
of the abundance difference A(//) - A shown in the bottom panels 
of each Figure. Here At//) is the value of A(O) obtained from fit¬ 
ting the line profile separately at the individual //-angle, while A 
denotes the A(O) value derived from simultaneously fitting the 
line profiles at all //-angles. 
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Fig.E.l. Same as in Fig. [3] except for the averaged 3D model atmosphere with <f m j c = l.Okm/s. 
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Pereira's spectra, Holweger-Mueller model atmosphere, E . =1.2 km/s 
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Fig.E.2. Same as in Fig. [3] except for the Holweger-Muller model atmosphere with E m j c = 1.2 km/s. 
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WCLC spectra, LHD model atmosphere, a MLT = 1.0, ^ mir = 1.2 km/s 
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Pereira's spectra, LHD model atmosphere, a MLT = 1.0, £ mic = 1.2 km/s 
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Fig.E.3. Same as in Fig. [3] except for the LHD model atmosphere with aMLT = 1-0 and £ m ; c = 1.2 km/s. 
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Fig.F.l. Same as in Fig. [7] except for the averaged 3D model atmosphere with 4 nnc = l.Okm/s. Best-fit solutions are not shown 
here, since they would require negative .S'h- 




Fig. F.2. Same as in Fig. |7] except for the Holweger-Miiller model atmosphere with 4mi c = 1.2 km/s. 
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Fig.F.3. Same as in Fig. [7] except for the LHD model atmosphere with (Tmlt =1.0 and <f ITIlc = 1.2 km/s. Best-fit solutions are not 
shown here, since they would require negative .S'n. 





























